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ABSTRACT 

We measure the location and evolutionary vectors of 69 Herschel -detected broad-line active galactic 
nuclei (BLAGNs) in the Mbh~ M* plane. BLAGNs are selected from the COSMOS and CDF-S fields, 
and span the redsliift range 0.2 < z < 2.1. Black-hole masses are calculated using archival spectroscopy 
and single-epoch virial mass estimators, and galaxy total stellar masses are calculated by fitting the 
spectral energy distribution (subtracting the BLAGN component). The mass-growth rates of both the 
black hole and galaxy are calculated using Chandra/XMM-Newton X-ray and Herschel far-infrared 
data, reliable measures of the BLAGN accretion and galaxy star formation rates, respectively. We 
use Monte Carlo simulations to account for biases in our sample, due to both selection limits and the 
steep slope of the massive end of the galaxy stellar-mass distribution. We find our sample is consistent 
with no evolution in the Mbh — M* relation from z ~ 2 to z ~ 0. BLAGNs and their host galaxies 
which lie off the black hole mass — galaxy total stellar mass relation tend to have evolutionary vectors 
anti-correlated with their mass ratios: that is, galaxies with over-massive (under-massive) black holes 
tend to have a low (high) ratio of the specific accretion rate to the specific star formation rate. We 
also use the measured growth rates to estimate the preferred AGN duty cycle for our galaxies to evolve 
most consistently onto the local Mbh — Mb u i relation. Under reasonable assumptions of exponentially 
declining star formation histories, the data suggest a non-evolving (no more than a factor of a few) 
BLAGN duty cycle among star-forming galaxies of ~ 10% (ler range of 1 — 42% at z < 1 and 2 — 34% 
at z > 1). 

Subject headings: cosmology: observations-galaxies: active-galaxies: evolution-quasars: emission lines 


1. INTRODUCTION 

As the most luminous persistent sources in the 
Universe, active galactic nuclei (AGNs), which are 
powered by the mass accretion of super-massive 
black holes (SMBHs), are widely believed to play 
important roles in the formation and evolution of 
typical massive galax ies (for a recent review, see 
iKormendv fe Hoi 12013ft . Indeed, over the past six¬ 
teen years, many studies have revealed that there are 


tight correlations between SA 
properties of host galaxies 

BHs and the physical 

(e.g., IMagorrian et al. 

1998 

: iFerrarese & Merrittl 

2000; Gebhardt et al. 

2000 

: Trema ne et al. 20021 

Marconi & Hunt 2003k 

Haring & Rix 2004; Giiltekin et al. 20091!. For example, 


the SMBH mass, Mbh, correlates well with the stellar 
mass of the bulge, Mb u i- The intrinsic scatter of the 
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Mbh — Mbui relation is found to be ~ 0.3 dex. Such a 
tight correlation might be explained by AGN feed back 
(e.g., iSilk fc Reed 119981 : iKind 120031 : iFabianl 12012ft . In 
this scenario, winds or jets launched from the accretion 
disk of the SMBH heat the interstellar medium (ISM) 
or clear the ISM out of the host galaxy. Star formation 
is therefore quenched due to the lack of cold gas. Such 
a coupled AGN triggering and feedback-regulated star 
formation could be responsible for the Mbh — Mbm 
relati on (see, e.g.. iDi Matteo et all [2005: Hopki ns et all 
120061! . On the other hand, a non-causal origin of 
the Mbh — Mr„i relati on has als o been pro posed by 
I Jahnke fc Macciol (120111! (see also I Pend 120071! . In this 
scenario, Mbh and Mb u i are not initially correlated. 
The Mbh — Mb u i relation is established simply by 
galaxy-galaxy mergers via the central limit theorem. 
Direct observational measurement of the evolution of 
the Mbh — Mb u i relation is crucial to investigate the 
formation and co-growth of SMBHs and their hosts. 
Observationally, the evolution of the Mrh — Mr„i re¬ 


lation is still a controversial topic (e . g..lPeng et al 
Treu et aT1l2007t I Jahnke et al.ll2009: Merloni et, al 


2Q06 


2010 ; 


Bennert et al. [20111) Schramm fc SilvermarJ2013ll . More- 

over, at high redshift (e.g., 2 ~ 1), it is very difficult to 
separate the bulge from the total stellar mass. Instead, 
most works explore the Mbh — M* relation, where M* is 
t he ga laxy total stellar mass. For instance, I Jahnke et all 
( 2009 ) find that at z ~ 1.4 the Mbh — M* relation is con¬ 
sistent with the local Mbh - Mb u i relation. Note th at the 
host galaxies of the local sample of lHaring fc Rbd (|2004f ) 
(hereafter HR04) are bulge dominated (i.e., Mb u i ~ M*). 
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Therefore, the I Jalinke et al.l (|2009T) result suggests little 
evolution of the Mbh — M* relation to z ~ 1.4. On the 
other hand, the hosts of the I Jalinke et al.l (l2009h sample 
have substantial disks (i.e., Mb u i < M*). That is, the 
Mbh ~ M Bu i relation does evolve as a function of red- 
shift. Making the l ocal Mrh — M b,,! relation then, as 
pointed out by iJahnke et all (120091 1, requires a redistri¬ 
bution of stellar mass from disks into bul ges i n AGN host- 
galaxies. Recently, the work of iSchramm fe Silvermaul 
( 2013 ) with a larger and less-bia s ed sa mple also sup¬ 
ported _Rie_£esults_of|jahnkeIer^I] ([2009) . On the other 
hand, iMerloni et alJ (|2010H compiled a broad emission- 
line AGN (BLAGh0) sample with well-measured multi¬ 
band photometry from the COSMOS survey. With 
this sample, they find that even the Mbh — M* rela¬ 
tion evolves (at 5c significance level) over cosmic time 
(i.e., the SMBH is over -massive at high redshifts; see 
also iBenner t et all 1201 lh . However, as pointed out by 
ILauer et al.l (1200^) . high-redslrift samples are generally 
selected by AGN activity (which depends on Mbh) and 
suffer from Eddington bias. The physical reasons are: (1) 
there is intrinsic scatter in the Mbh — Mbh or Mbh — M* 
relation; (2) the galaxy stellar mass function is bottom 
heavy (i.e., the number density of galaxies steeply de¬ 
creases to higher stellar mass). Hence, for these high- 
redshift samples, an AGN with given Mbh is more likely 
to be found in less massive galaxies, resulting in an ap¬ 
parent evolution of the Mbh — Mb u i or Mbh — A/* re¬ 
lation (also see Section 14.11) . This bias increases with 
redshift as at high redshifts we often only have access 
to very luminous AGNs (very massive SMBHs). Addi¬ 
tionally, Mbh estimation itself is also biased for these 
high-redshift samples (|Shen fc Kellvll20Tol ). The reasons 
are similar: (1) Mbh is estimated from AGN luminosity 
which also has intrinsic scatter; (2) the AGN luminos¬ 
ity functio n is also bottom he avy. Note that the sam¬ 
ple of ISchramm fc Silverman (120131 ) consists of mostly 
l ower luminosity AGNs and hence suffers less bias than 
IMerloni et al. ( 201 (f). Bias corrections must be made to 
compare the high-redshift Mbh — M* relation with the 
local one (for a more recent di scussion on this topic, see, 
e.g., iSchulze fe Wisotzkill2014D . 

In addition to measuring the Mbh — M* relation at 
high redshift, it is also valuable to study its instanta¬ 
neous evolution; i.e., the SMBH and galaxy stellar mass 
growth rates. Especially, will future SMBH and galaxy 
mass increases help establish/maintain the Mbh — M* 
relation? This question has been indirectly investigated. 
For example, many works explored the possible connec¬ 
tion between star formation rates (SFRs) and AGN lumi- 


nosities (e.g.. Raffertv et al. 

2011: Mullanev et al.l 20121 

Rosario et al. 2012, 2013aib 

; Chen et al.l 20131 and also 


see Alexander & Hickox 2 012 for a re view). However, 
as pointed out by iHickox et al.l (|2014f) . AGN variabil¬ 
ity can dilute the intrinsic correlation between star for¬ 
mation and AGN activity. Therefore, this connection is 
expected to be observed only in an average sense, as in¬ 
dividual galaxies will tend to show a much weaker (or 
absent) correlation. The average SMBH accretion rate 
traces the average SFR so well that the local Mbh — M* 
relation can be established. 

8 Throughout this work, we use “BLAGN” and “AGN” inter- 
changeably. 


Simultaneously exploring both the Mbh — M* rela¬ 
tion at high redshifts and its instantan eous ev olu tion 
is eve n more valuable but is challenging. IMerloni et al.l 
( 2010 ) explore this issue by measuring Mbh, M*, and 
their growth rates. They find that star formation and 
AGN activity can actually help reduce the scatter of the 
Mbh — M* relation. However, the SFR in their work is 
estimated by fitting a stellar population synthesis model 
to the optical-to-NIR spectral energy distribution (SED) 
which has an uncertainty of ~ 0.7 dex. A more accurate 
determination of SFR (e.g., with Herschel observations) 
could better elucidate the coupled growth of AGNs and 
their host galaxies. 

In this work, we uniquely explore the Mbh — A/* re¬ 
lation and its evolution, accounting for biases and using 
robust and uncontaminated estimates of SMBH growth 
and galaxy star formation rate. In Section [2j we de¬ 
scribe our sample in detail. In Section^ we describe our 
methods for estimating Mbh, A/*, M, and SFR. In Sec¬ 
tion |4] we discuss selection biases and present our results. 
In Section EH we discuss implications of our results. Fi¬ 
nally, we summarize our main results and outline the 
possible future improvements in Section 0 Through¬ 
out this work we adopt a flat ACDM cosmology with 
H 0 = 70 km s -1 Mpc -1 , 9 m = 0.3 and 9 a = 0.7. 

2. SAMPLE CONSTRUCTION 

The 2-deg 2 Cosmic Evolution Survey (COSMOS, 
iScoville et alJ I2007D and 464.5 -a rcmin 2 Chandra Deep 
Field-South surv ey (CDF-S; e.g.. iLuo et al.ll2008l 120101 : 
iXue et ahl 1201 lh both have well-calibrated multi-band 
photometry (from the X-ray to radio bands) and high- 
quality spectroscopic coverage. Hence, we selected our 
sample from the COSMOS and CDF-S fields. First, we 
selected BLAGNs based on X-ray detection of AGN ac¬ 
tivity (f or t h e advantages of using X-ray selection, see 
iBrandt fc Alexander! 12010 ) and subsequent optical spec¬ 
troscopic identification of broad emission lines (with cov¬ 
erage of either H/3 or Mgn). Second, we cross-matched 
our BLAGNs with the Herschel PEF0 dLutz et al.ll201 lh 
(for the CD F-S field, we also used the GOODS -Her schel 
catalogs, se elElbaz et al.ll20llt lMagn elli_et_a J1l2013l) and 
HerMEf£3 (lOliver et al.l l2012t iWang et~ahl 120131 ) cata¬ 
logs. Third, the new sample was cross-matched with ex¬ 
isting multi-band photometric catalogs to build the final 
sample. In the following subsections, we discuss our sam¬ 
ple selection in detail (see also Table 1). 

2.1. COSMOS Field BLAGNs 

The 2-deg 2 COSMOS field has been surveyed in 
the X-ray by XMM-Newton dCappelluti et all 120091 : 
iBrusa et al.l l2010h and partialh^ (0.9 de g 2 ) by Chan¬ 
dra (lElvis et al.ll2009l : iCivano et, al.ll2012D . We first se¬ 
lected all z < 2.4 BLAGNs (ensuring the coverage of ei¬ 
ther H/3 or Mg il) f rom the public COS MOS spectroscopy 
dLillv et al.l l2007t I Trump et al.l l2009afl . Our BLAGN s 
have public Magellan/ IMACS dTrump et al.l l2009al) . 
SDSS (I York et alJl2000l) . or VLT/VIMOS (zCOSMOS; 

9 http://www.mpe.mpg.de/ir/Research/PEP/index 

10 http://hedam.lam.fr/HerMES/ 

11 This together with the on-going Ch andra COSM OS Legacy 
Survey will cover the full COSMOS field (jCivano et al.l[2013n . 
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iLillv et aLll2007f) spectroscopic datrQ which allow us to 
estimat e M bh via the single-epoch approach (see also 
Section 13.11) . The Ma gellan/IMACS spectr oscopy is lim¬ 
ited by «ab < 22.5 (jTrump et, al.l l2009al) . For consis¬ 
tency and ease of modeling the observational biases, we 
impose the same iab < 22.5 limit for the combined COS¬ 
MOS and CDF-S sample. We then cross matched them 
with X-ray catalogs to get the corresponding (rest-frame) 
2 — 10 keV luminosity. For X-ray observations, we pre¬ 
fer Chandra over XMM-Newton and hard band over soft 
band (i.e., the Chandra hard band has the highest prior¬ 
ity, and the XMM-Newt on soft band has the lowest pri¬ 
ority). As pointed out bv lTrump et al.l (l2009al l20TTTi and 
will be shown in Figure [3] our AGN selection is limited 
by the /ab <22.5 criteria and not the X-ray sensitivity 
(whether XMM-Newton or Chandra). 

To get a robust estimate of SFR, we cross-matched 
our BLAGN sample with the Herschel PEP (100 /im 
and 160 /tm) and HerMES (250 /mi) catalogs. Note 
that for the HerMES catalog the SPIRE/ Herschel ob- 
ser vat ions have a large point-spread function. We fol¬ 
lowed iChen et al.1 12 0131 ) and chose a maximum match¬ 
ing radius of 5 . The Herschel PEP catalogs have im¬ 
proved astrometry because sources were extracted using 
Spiteer/MIPS 24 /mi priors, and so we used a maximum 


matching radius of 2 '0 As for multi-band photometry, 
we us ed the COSMOS P hotometric Redshift Catalog Fall 
2008 Albert et al. 200 0; to get the optical-to-NIR broad¬ 
band photometry information ( B , V, g , r, 2 of Subaru; u, 
i, Ns of CFHT; and J of UKIRT) for the selected so urces 
(see also iCapak et al .1120071 iMcCracken et al.ll2010D. We 
also i ncluded four Spitzer/ IRAC channels (j Sanders et al.1 
120071 ) and GALEX /Near UV d ata when ava i lable. Zero- 
point corrections suggested in lllbert et~. (feOOfl ) were 
applied. The Galactic extinction was also corrected us - 
ing the Galactic extinction map of lSchlegel et al.l (|1998t ). 
Note that for the optical-to-NIR bands, the Galactic ex¬ 
tinction corrections are already included in the catalog. 


2.2. CDF-S Field BLAGNs 

The 464.5-arcmin 2 CDF-S field is covered by the 
deep est X-ray survey, with 4 Ms of Chandra cover¬ 
age dXue et al.1 1201 ll) . We selected BLAGNs by cross¬ 
matching t he 4 Ms CDF-S catalog w ith the optical spec¬ 
troscopy of iSzokolv et al.l (|2OO4D0 Similar to the COS¬ 
MOS sources, we only selected BLAGNs (with »ab < 
22.5) whose host galaxies are also detected by the Her¬ 
schel PEP or HerMES (the cross-matching criteria is the 
same as that of COSMOS BLAGNs). As for multi-band 
photometry, we used the MUSYC0 32-band (including 
U, U 38, B , V, R, /, 2 , J, H , K bands; 18 Subaru 
me dium bands; and four S pizter /IRAC channels) cata¬ 
log dCardamone et aLll2010t and references therein). The 
G alactic extinction and zer o-point corrections suggested 
in iCardamone et al.l 12010 1 were adopted. 


12 http://irsa.ipac.caltech.edu/Missions/cosmos.html 

13 For the Herschel catalogs, our maximum matching radii lead 
to a false match rate of 1.3% for PEP and 9% for HerMES. (All but 
five objects have both PEP and HerMES data.) The false-match 
ra te is calculated by following the procedure described in Section 3 
of ll.no el all 1)2008 1. 

14 http: //tow. mpe.mpg.de/CDFS/data/ 

15 http: //tow. astro.yale.edu/MUSYC/ 


2.3. Total BLAGNs 

Based on the above criteria, we started with a total 
number of 251 X-ray BLAGNs at 2 < 2.4 (ensuring cov¬ 
erage of either H/3 or Mgil), with 221 from the COSMOS 
field and 30 from the CDF-S field. When cross-matched 
with the Herschel catalogs, 96 sources were left (85 from 
the COSMOS field; 11 from the CDF-S field). 69 of 
these 96 sources have high-quality spectroscopic obser¬ 
vations and robust FIR detections which enable us to 
measure Mbh, M*, M, and SFR, reliably. For the other 
27 sources, either the spectroscopic data (nine objects) 
or the FIR observations (18 objects) are not sufficient for 
a relia ble measurement of A7 bh or SFR (see Sections 13.11 
fc 13.41) . Note that due to the Herschel sensitivity, ~ 2/3 
of BLAGNs are undetected. We will consider the bias 
introduced by the Herschel sensitivity in Section 14.11 

For the 69 sources in our sample, 62 are from the COS¬ 
MOS field and seven are from the CDF-S field. Note that 
their spectra are flux calibrated. We compared the phys¬ 
ical properties (Section [3| of these seven CDF-S sources 
with the COSMOS sources (using the Mann-Whitney U 
test) and found no statistical difference (i.e., the null 
probability p > 0.05). 

3. MEASURING PROPERTIES OF AGNS AND 
THEIR HOST GALAXIES 

3.1. Black-Hole Mass Estimation 

We adopted single-epoch virial SMBH mass estima¬ 
tors to estimate SMBH masses. More specifically, H/3 
or Mg ll line estimators were used depending on the red- 
shift. These SMBH mass estimators use the correlation 
between the radius of the broad line region (BLR), Rblr, 
and the continuum luminosity, A L\, Rblr oc (ALa)° , 
reveale d by reverberation mapping observa tions of local 
AGNs dBentz et al.ll2006t iKasni et al.ll2007t ). If the BLR 
is virialized, then the SMBH mass can be estimated by 
M B h = /Rblrt|whm> where dfwhm is the full width 
at half maximum of the broad emission line and / is re¬ 
lated to the BLR geometry (calib rated from dynamical 
estimates of Mbh, see, e.g.. lOnken et al.ll2007t ). There¬ 
fore, we can estimate the SMBH mass from A L\, and 
ffwhm, 

j|/fnTT 

log ( WT - ) = A + B log(AL A ) + 2 log(uFWHM), (1) 
Mq 

where A L\ and cfwhm are in units of 10 44 erg s~ 4 
and 1000 km s _1 , respectively. This rela tion has an 
intrinsic scatter of ~ 0.4 dex dVestergaard fc Petersonl 
I2006D . which is due to the unknown BLR geometry and 
other factors (e.g., ionization state). Note that for H/3, 
A = 5100 A, A = 6.91, B = 0.5 dVestergaard fe Osmerl 
l2009i): for Mg II, A = 3000_A, A = 6.86, B = 0.47 
dVestergaard fc Petersonl I2006D . Both A L\ and ufwhm 
can be measured by fitting the single epoch spectra. 

We performed an iterative chi-squared nrinimizatior0 
to fit the broad emission lines (H/3 or Mgil) for each 
source and measured vfwhm and A L\. The fit t ing pro - 
cedures were similar to those of I Trump et al.1 (l2009bD . 
That is, a pseudo-continuum (power-law continuum plus 

16 Using the Imfit python package, available from 
http://cars9.uchicago.edu/software/python/lmfit/index.html 
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broadened Fe template of IVestergaard fe Wilkesl 12001(1 
and one or two broad Gaussian line profiles were (simul¬ 
taneously) used to get a best fit of each spectrum. For 
H/3, we also added [O ill] AA4959, 5007 lines and a narrow 
H/3 component to the fit and then subtracted them when 
fitting the broad emission lines since these narrow com¬ 
ponents come from the narrow-line region. Our criterion 
for the use of one or two broad Gaussian line profiles was 
as follows: first, we used two Gaussian components to fit 
the broad emission line and then we compared ufwhm 
of both Gaussian components. If at least one component 
had vfwhm < 10 3 km s , we refitted the spectrum with 
only one Gaussian component. Note that if we use two 
Gaussian line profiles for all sources, our Mbh estimates 
change by no more than 0.2 dex, and this is small com¬ 
pared to the intrinsic uncertainty of Mbh- 

In Figure CD we present examples of our fits to the 
broad emission lines. The gray and red solid lines rep¬ 
resent the observed spectra and our best fit, respec¬ 
tively. Red and orange dot-dashed lines represent the 
power-law continuum component and Fe template, re¬ 
spectively. The dotted lines correspond to minor fea¬ 
tures (i.e., [O ill] AA4959,5007 lines, the narrow H/3 line, 
which are excluded from the broad-line fit). Blue dashed 
lines represent best fits of broad emission lines. We cal¬ 
culated vfwhm as follows: if two broad Gaussian profiles 
were used, we calculated ufwhm from the emission pro¬ 
file which is the summation of the two broad Gaussian 
components; otherwise, we calculated vfwhm from only 
the one broad Gaussian profile. We have visually in¬ 
spected all fitting results and found nine sources whose 
spectroscopic line profiles do not allow reliable measure¬ 
ment of Mbh o We excluded these nine sources from our 
sample. Meanwhile, A L\ was calculated from the power- 
law continuum component. Note that we also applied a 
small correction to account for the contamination of the 
stellar light of the host galaxy suggested in Section 13.21 

A summary of SMBH masses can be found in Ta¬ 
ble 2 (before compiling the Table, we removed another 
18 s ourc es because of inadequate FIR data, see Sec¬ 
tion [T4|)- There are 24 (46) sources in Table 2 whose 
Mbh were m easured using H/3 (Mgn). Previous stud¬ 
ies (e.g., iShen fe Liul l2012p have revealed that the two 
A7 bh estimators are consistent. Four sources in Table 
2 have spectroscopic data enabling us to measure Mbh 
using both H/3 and Mgn. For three BLAGNs, the two 
Mbh estimators are consistent with each other within 0.2 
dex. For the remaining BLAGN (ID-48), the two Mbh 
estimators show a larger deviation (0.6 dex) but are still 
consistent within the 2er uncertainty. In any case, we 
preferred the H/3 estimator over Mgn estimator. 

3.2. Galaxy Total Stellar Mass Estimation 

We utilized an SED fitting technique to estimate M*. 
Since we were dealing with BLAGNs, the SEDs (from 
UV to Near IR) of the host galaxies are contaminated 

17 Four of these nine sources only show narrow H/3 lines in their 
spectra; i.e., they are misclassified as BLAGNs, and are instead 
Type 2 or intermediate type (type 1.x) AGNs. Five of these nine 
sources show strong absorption features superposed on the broad 
emission lines. The removal of these sources should not introduce 
a selection bias to our sample, since they are essentially random 
occurrences which are unlikely to be strongly correlated with Ed¬ 
dington ratio or host properties. 


by AGN emission. We th erefore fo llowed the approach 
of iBongiorno et all (|2012T ) (see also lMerloni et al1l2010l) 
which fits the observed SEDs as a composite of both AGN 
and galaxy stellar emission. That is, 

/obs,A = CaGn/aGNA + CGal/Gal.A) (2) 

where /agn.a and /Gai.A are the fluxes from the AGNs 
and host galaxies, respectively. For SEDs, fourteen 
bands (observed-frame wavelength ranging from 2500 A 
to 79595 A) are used for the COSMOS sources and 32 
bands (observed-frame wavelength ranging from 3656 A 
to 79595 A) are used for the CDF-S sources. For 
A GN emissio n, we adop ted the mean SED template 
of iRichards et all ( 2006 ). For galaxy templates, we 
selected from an S EP l ibrary constructed using the 
iBruzual fe Chariot] (120030 stellar population synthesis 
model. We ass umed an initial mass function (IMF) of 
[Chabrierl (120031 ) and adopted ten exponentially declining 
star formation histories (SFHs), i.e., SFR oc e _tage,/T , 
with star formation timescale r ranges from 0.1 to 
30 Gyr, plus a SFH with constant SFR. The galaxy age 
ranges from 50 Myr to 9 Gyr with an additional con¬ 
straint that the galaxy age should not be larger than 
the age of the Universe at the redshift of the source. 
The SED library was generated iterating over r and the 
galaxy age. When fitting the observed SED, we also took 
intrinsic extinction into consideration. For AGNs, we 
used a Small Magel lanic Cloud-like dust-reddening curve 
(|Prevot et al.l [l984) with E(B — V) < 1.0 For galaxies, 
we ad opted the Calzetti extinction curve dCalzetti et al.1 
12000(1 with E(B — V ) < 0.5 if t^/r < 4, otherwise 
E(B — V) < 0.15 (iFontana et al .1120061 : iPozzetti et al.1 
l2007t IBongiorno et al.ll2012fb 

We additionally required that the AGN continuum 
emission should be larger than 10% of the galaxy contin¬ 
uum emission at 2800 A (if the Mgn virial mass estima¬ 
tor was used) or 4861 A (if the H/3 virial mass estimator 
was used). This was simply motivated by the fact that 
our sources all have observed broad emission lines (H/3 or 
Mgil) indicative of a significant AGN contribution. We 
chose 10% based on the following considerations: (1) we 
performed a simple simulation by creating a mock spec¬ 
trum with a young galaxy SED, an AGN (power-law) 
component, a broad Gaussian emission line with a typical 
equivalent w idth of 30 A (60 A) for Mgil (H/3) (e.g., Sec¬ 
tion 1.3.4 of lPetersonlll997D . and white noise (assuming 
a signal-to-noise ratio of ~ 10); (2) we refitted the mock 
spectrum with a power-law and a Gaussian emission line 
and found that we cannot reliably fit the emission line 
if the AGN emission is less than ~ 10% of galaxy emis¬ 
sion. Furthermore, we have performed a Monte Carlo 
sim ulati on to mimic our selection procedures (see Sec¬ 
tion [4J] for more details). For our simulated sample, we 
also found that there are a negligible number of sources 
with AGN emission less than 10% of the galaxy emission 
at 2800 A (or 4861 A). 

We fitted f 0 bs,\ to the observed SEDs by perform¬ 
ing iterative (reduced) chi-squared minimization. Fig¬ 
ure [2] shows examples of our two-component (AGN plus 
galaxy) fitting to the observed SEDs. The data are well 
described by a combination of AGN and galaxy emis¬ 
sion. We also inspected the corresponding HST/ACS 
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images and found that generally our SED fits are consis¬ 
tent with the morphology information. That is, when the 
SED fit suggests the AGN emission dominates in the i 
band, the HST/ACS image also shows a point-like mor¬ 
phology (the inverse is also true). From each best fit, 
we determined both the normalization, t age , r, and the 
SFH of our galaxy SEP . We then used them and the 
iBruzual &: Chariot! (I2003f) model to calculate the galaxy 
total stellar mass. A summary of the galaxy total stel¬ 
lar masses can be found in Table 2. Note that the two- 
component SED fit can also tell us the AGN contribution 
to either the 3000 A or 5100 A luminosity, and we can use 
this information to calculate the true AGN contribution 
to A L\ measured in Section EU This small correction 
(accounting for the contamination of the stellar light of 
the host galaxy) has been included when calculating A L\ 
from the continuum component in Section 13.11 


3.3. Accretion-Rate Estimation 

To measure the accretion rate, M, we esti mated the 
bolometric luminosity, TboI fe.g.. lSoltanlll982D . The ac¬ 
cretion rate is, 


M = (1 ~ v)Lb ° 1 

nr^ 


( 3 ) 


where ry = 0.1 is the assumed radiative efficiency of the 
accretion disk, and c is the speed of light. We used two 
methods to calculate TboU (1) the rest-frame 2—10 keV 
luminosity alone; (2) the SED fit plus the rest-frame 2 — 
10 keV luminosity. 

We first used the rest-frame 2—10 keV luminosity as 
an estimator of LboI- The rest-frame 2—10 keV luminos¬ 
ity is calculated from the observed-frame 2 — 10 keV flux 
(assuming a typical power-law photon index of T = 1.7). 
Note that there are four sources in our sample that are 
undetected in the observed-frame 2 — 10 keV band. For 
them, we instead calculated the rest-frame 0.5 — 2 keV lu¬ 
minosity from the observed-frame 0.5 — 2 keV flux (again 
with T = 1.7) and used (different) bolometric corr e ctions 
to estimate LboI- We adopted the lHopkins et al.1 (|2007f l 
luminosity-dependent bolometric corrections: 


^Band 


L Bol _ , , TboI , fcl , / TboI w - 2 
- -“Hminr ) ' “ 2 1775777 ) ’ 


L Band 


10 10 T© y 


'lO 1O I/ 0 ' 


( 4 ) 


where «Band is the bolometric correction fac¬ 
tor. The constants (di, k\, d%, fe) are given 

by (17.87,0.28,10.03,-0.02) for To. 5-2 keVi and 
(10.83,0.28, 6.08,—0.02) for £ 2- 10 keV- The uncertainty 
of KBand is (Hop kins et al.ll2007li . 


01og( K ) = CTATBol/TO 9 ^©)^ + (72, (5) 


with (ex!, /3, (j 2 ) = (0.046,0.10,0.08) in the 0.5 — 2 keV 
band and = (0.06, 0.10, 0.08) in the 2 — 10 keV band. For 
sources in our sample, the median of the uncertainties is 
~ 0.2 dex. 

We also combined the SED fitting results presented 
in Section 13.21 with X-ray observations to calculate the 
bolometric luminosity. This pr ocedure was similar (but 
not identical) to that of lTrump et al.l (|2011[ ). First, we 
used the rest-frame 2—10 keV luminosity to get the 
0.5 — 250 keV luminosity (assuming T = 1.7 and a 
cut-off at 250 keV); second, we integrated the big blue 


bump component (from 30 A to 10 4 A) in the best-fitted 
AGN SED which should be responsible for the emission 
from the accretion disk (note that we neglected the IR 
“torus” emission as it is probably reprocessed); finally, 
we obtained the bolometric luminosity by summing the 
0.5 — 250 keV luminosity and the total big blue bump lu¬ 
minosity. This approach gives a direct estimate of TboI- 

We compared the two bolometric luminosity estimates 
and found that the ratios of these two estimates follow 
a log-normal distribution. This log-normal distribution 
has a mean value of 0.03 and a standard deviation of 
0.4 dex. This deviation should be caused by the un¬ 
certainties of the two bolometric luminosity estimators. 
Therefore, we conclude that the uncertainty of the bolo¬ 
metric luminosity estimated using the SED fit plus the 
X-ray luminosity is (0.4 2 — u^,^) 0 ' 5 = 0.35 dex. In the 
following calculations, we use the bolometric luminosity 
calculated from SED fits plus X-ray observations to mea¬ 
sure M. A summary of the bolometric luminosities can 
be found in Table 2. 

Figure [3] plots LboI as a function of redshift. In ad¬ 
dition, we also include the bolometric luminosity limits 
introduced by the *ab < 22.5 spectroscopy limit and 
the X-ray sensitivity limits of COSMOS (soft bands of 
Chandra and XMM-Newton) and CDF-S (hard band of 
Chandra ). The X-ray limits are significantly deeper than 
the optical spectroscopy limit, and the sample is essen¬ 
tially limited only by iab < 22.5 (see also iTrump et al.l 
l2009al) . The four sources which are only detected in 
the observed-frame 0.5 — 2 keV band are denoted as red 
points. The redshifts of our BLAGNs span 0.2 < z < 2.1 
with a median redshift of 1.18. We also define a “high- 
redshift” sub-sample of the 45 sources at z > 1, and a 
“low-redshift” sub-sample of the 25 sources at z < 1. 

Figure [4] plots the bolometric luminosity as a func¬ 
tion of M B h- We also include the lines of LboI = 
0.01 TEdd, Tsoi = 0.1 TEdd and LboI = ^Edd, where 
Tfidd = 1-26 x 10 38 M B h/A4© erg s _1 . As clearly seen 
from this figure, the Eddington ratios of our BLAGNs 
are typically spread from 0.01 — 1 (consistent with, e.g., 
Kolhneier et all 120061 : ITrump et al.l I2009IJ : ILusso etall 
20121) . For COSMOS sources, we also made a detailed 
compari son between Edding ton ratios of our sources and 
those of lTrump et al.l (1201111 by checking the overlapping 
sources one by one. We found that the median absolute 
deviation (M.A.D.) of the differences is 0.2 dex, much 
smaller than the uncertainty we ass igne d to Mbh (see 
Section 13.51) . For the CDF-S field, iBabic et all (1200711 
found a broad distribution of Eddington ratios, using 
host galaxy mass (or velocity dispersion) to estimate 
SMBH mass. Of their sources with broad emission lines, 
6/11 have Eddington ratios between 0.01 and 1, and 
two others have Eddington ratios consistent with 0.01. 
Only 3/11 have Eddington ratios significantly smaller 
than 0.01. Despite the very different methods of Mbh 
estimation, our results are broadly consistent with their 
results. 


3.4. SFR Estimation 


Far-IR emission has long been argued to be an excellent 
SFR tracer as it is larg ely free of dust e xtinction. We 
used the widely adopted iKennicuttl (|1998) relation (with 
a modification to account for the IMF of lChabrieill2003H 
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to derive the SFR: 

SFR(M 0 yr _1 ) = 1.09 x 10 _10 Air/A q , (6) 

where Air is the total 8 — 1000 pm IR luminosity. 
BLAGNs make significant contributions in the near- to 
rnid-IR bands. The FIR band is, how ever, known to be 
domi nated by galaxy emission (e.g., [Kirkpatric k et al.l 
120121 hereafter K12). In this work, we used FIR data 
from the Herschel PACS/SPIRE bands to normalize 
the z ~ 1 star-forming galaxy SED template of K12 
and integrated the normalized SED to estimate Fir. 
There are other popular galaxy IR SEP templates (e.g., 
iCharv fe ElbazlDOOll iDale fe Heloull2002D . However, the 
K12 template is directly derived from the high-redshift 
sources and therefore may be more relevant to our high- 
redshift sources. As stated before, we have additionally 
removed 18 sources which have z > 1 but are detected 
only at the 100 /im band. The reason is that, under such 
circumstances, the 100 pm band might be sig nificant ly 
conta minated by AGN emission (e.g., K12. lNordon et ahl 
120121 1. Such measurements of Air are therefore not re¬ 
liable. For our 69 sources, 26 are detected at 100 pm, 
160 /im, and 250 pm. 24 are detected in two bands. The 
remaining 19 sources are detected only in one band. A 
summary of Air can also be found in Table 2. 

To check the consistency of our Air estimation, we 
took our 26 sources that are detected in three bands 
and compared Air values that are estimated from all 
three bands with those values from two bands or one 
band. Figure [5] shows the result. As seen from the fig¬ 
ure, Air values estimated from three bands are consistent 
with those from only two bands or one band. There¬ 
fore, we can conclude that our method for estimating 
Air is self-consistent. Note that, even for 2 > 1 sources, 
Air estimated from 100 /im agrees well with that from 
three bands, which indicates that for these three-band 
detected sources AGN emission does not strongly con¬ 
taminate even the blue Herschel FIR band. For the re¬ 
moved 18 sources, since we only have information on the 
bluest band (i.e., 100 (im) which corresponds to rest- 
frame 30 — 50 /im (for z = 1 ~ 2), we cannot rule out 
the possibility that AGN contamination may still be im¬ 
portant [3 

To quantify the star formation activity of our sam¬ 
ple, we compared the SFR — M* relation of our galaxies 
with t he star-forming “ma i n sequence” (e.g . .lElbaz et al.1 
120071 iNoeske et al.1 120071 : I Whitaker et al.1I2012D . The 
SFR — M* relation (i.e., t he star-forming “main se¬ 
quence” ) we adopted is from IWliitaker et al.1 (1201211 , 

log SFR = 0 ( 2 ) (log M* - 10.5) + b(z), (7) 

with a(z) = 0.7 — 0.13,3, b(z) = 0.38 + 1.143 — 0.192 2 . 
This relation has a scatter of 0.34 dex, roughly indepen¬ 
dent of z and M*. Figure [6] plots the ratio of the SFR 
we measured to the one predicted by Equation [7] For 
the uncertainty of the ratio, we combined the uncertain¬ 
ties of M*, SFR, and the scatter of Equation [71 Our 
sources show higher star formation activity than that of 
the main sequence (although they are not significantly 
offset by more than their uncertainties; also the mean 
ratio is ~ 2, unlike starburst galaxies whose ratios have 

18 Note that our results below do not change materially if we 
include the 18 objects by assuming their SFR estimation is robust. 


> 3cr offset from Equation 0 . This is driven mostly by 
the Herschel flux limit, which lies aroun d or abov e the 
main sequence (see, e.g., Figure 4 of lNordon e t al] l2012f) 
and causes an Eddington bias (see Section [XT]). The fac¬ 
tor of two excess in the mean ratio is fully consistent 
with the expected Eddington bias if the host galaxies 
of our BLAGNs were drawn from the star-forming main 
sequence. 


3.5. Error Budget 


In our sample, Mrh, M*, M, and SFR all have signif¬ 
icant uncertainties. Let us first consider the uncertainty 
of Mbh estimation. For a fraction of our sources, Mbh 
has been carefully estimated by some previou s studies 
(e.g.. iMerloni et ahll201d : iMatsuoka et al.l[20l3l) . There¬ 
fore, we compared our Mbh estimations with these pre¬ 
vious works. We found that there is no global system¬ 
atic offset and the standard deviation is no more t han 
0.3 d ex (actually, when comparing with Matsuok a et al.1 
120131 we found that for all sources but one, the deviations 
are smaller than 0.2 dex; for the remaining one, the devi¬ 
ation is ~ 0.3 dex). Therefore, the error is dominated by 
the intrinsic error of the SMBH mass single-epoch virial 
relation. We then a dopted 0.4 dex, the intrinsic error, 
as the error of Mbh dVestergaard fe Petersonll2006l) . For 


M, we only considered the uncertainty from Ab 0 i, which 
is ~ 0.35 dex (see Section 

For Ms, we compared our results with those 
of IBongiorno et all (|2012H (COSMOS sources) and 
IXue et all (1201011 (CDF -S sources). We found the me¬ 
dian and M.A.D. of the differences is 0.02 dex and 0.2 
dex, respectively. Note that there are a handful of 
sources in our sample showing r elatively large dev i ations 
of M* with respect to those of IBongiorno et ah ( 201 2; 
(although there is no global systematic offset). This is 
likely caused by several factors, e.g., different photomet¬ 
ric d ata (our data included t he NUV data from GALEX 
while Bongior no et al.ll2012t on the other hand, took the 
24 pin data into consideration); our a dditi ona l re quire- 
ment of AGN contribution. Note that IBongiorno et al 


(2012) and our work both used the IMF of [Chabrier 


(2003), and the galaxy stellar masses from IXue et al 


1 2010 ) we re als o converted to those with the IMF of 
iChabried (|200l) . The M.A.D. we obtained should indi¬ 
cate this systematic/methodology uncertainty when esti¬ 
mating M* from SEDs. This uncertainty is much larger 
than the uncertainty introduced by the photometry error. 
Therefore, we adopted the normalized M.A.D. (NMAD), 
onmad = 1-5 x M.A.D. = 0.3 de x, as the uncertain ty of 
the galaxy total stellar mass (jMaronna et al.ir2006D . 

SFR estimation is subject to two major uncertainties: 
(1) the errors of the FIR fluxes; (2) the intrinsic uncer¬ 
tainty of using the K12 template to estimate SFR is, as 
reported by K12, 0.17 dex. We added these two uncer¬ 
tainties in quadrature to account for the full uncertainties 
of SFRs. 


4. CO-EVOLUTION OF AGNS AND THEIR HOST 
GALAXIES 


19 For simplicity, in this work, we n eglec t the scatter of tj. 
Note that our anti-correlations in Section 16.21 do not significantly 
change if we instea d assume the scatter of rj is ~ 0.5 dex (e.g., 
IDavis Hr. Laoiir201lh 
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4.1. Selection Biases 

We now discuss biases due to selection effects. First, 
our sources are selected based on AGN luminosity (as 
well as the luminosity contrast of the AGN and host 
galaxy in the optical/UV) and broad emission lines (H/3 
or Mgll). Selecting galaxies around luminosity-limited 
AGNs leads to an Eddin gton bias in_the Mbh — M * re¬ 
lation, as detailed by lLauer et al.l ([200 7]) (a lso see Sec¬ 
tion [T]). On the top of the lLauer et al-T i 20071 ) bias, there 
is anoth er bias on Mbh for our luminosity-limited AGN 
sample dShen fe Kellvl l20lol) . Second, we required our 
sources be detected by Herschel. The Herschel detection 
limit effectively introdu ces a bias which acts in the oppo¬ 
site direction as the lLauer et al i (j2007f ) bias, since a SFR 
limit generally selects more massive host galaxies for an 
AGN with given Mbh- 


4.1.1. The Basic Model to Estimate Biases 

To quantify selection biases in our sample, we be¬ 
gin by assuming that the scatter in each quantity has 
a lo g-normal distribut i on (as commonly a ssumed in, 
e.g.. lLauer et al.112007t IShen fe Kellvl 120101) . Accord¬ 
ing to the local SMBH mass-galaxy stellar mass rela¬ 
tion (to) = cis + C 2 , the distribution of the SMBH 
mass, Mbh = 10 m Mg, at fixed galaxy stellar mass, 
M* = 10 s M 0 , is 

p{m\s) = (27rcr M ) exp(--), (8) 

where is the intrinsic scatter. We also assume that our 
BLAGN host galaxies are drawn from the star-forming 
main sequence in which SFR correlates well with galaxy 
stellar mass. The motivations are as follows. First, we 
demonstrate below that this is a good assumption, as 
the apparent factor of two excess of SFR seen in Sec- 
tion !3.4l and Figure[G]is fully consistent with the b ias f rom 
the He rschel sensitivity limit. Moreover, iRosario et al.1 
(l2013bD also demonstrated that, after using simulations 
to model selection biases carefully, BLAGN hosts are 
consistent with being drawn from the star-forming main 
sequence relation. Note that BLAGNs are also observed 
to be a factor of > 5 less common in q uiescent galax¬ 
ies than in blue star- forming hosts dTrump et al.l 120131 : 
iMatsuoka et al.l 120141) . Therefore, we can connect SFR 
with galaxy stellar mass by Equation!?] The distribution 
of SFR (SFR = 10 A M 0 yr -1 ) for fixed s is 


g(\\s) = {2Ttcr 2 F ) 0 5 exp( 


(A — (a(s — 10.5) + b)) 2 


( 9 ) 

where a = 0.7 — 0.132:, b = 0.38 + 1.142 — 0.192 2 (see 
Equation [?| ■ The intrinsic scatter of the SFR — M* re¬ 
lation o\f is roughly ind ependent of both red shift and 
galaxy stellar mass ('e.g.. IWhitaker et al.ll2012D . 

Equations [8] and M can be used to estimate the bias 
in the Mbh — M* relation resulting from sample selec¬ 
tions based on Mbh and SFR. While the SFR selec¬ 
tion is directly related to the Herschel sensitivity, the 
Mbh selection is a more complicated function of the 
AGN luminosity limit and the broad-line width limit 
(i.e., ufwhm > 10 3 km s -1 ). M otivated by Shen et all 
( 2008). [Trump et all (1201 ll) . and lLusso et al. (|2012f) we 


related Mbh and LboI by assuming our BLAGNs have 
a log-normal distribution of the Eddington ratio cen¬ 
tered on 0.1 with a scatter of <tl dex and a cut off of 
< 0.01 and > 1. Using a different shape for the distribu¬ 
tion (e.g., log-uniform) turns out to have little effect on 
our bias estimate. However, assuming a different mean 
Eddington ratio for the distribution does have a small 
effect; see Section 14.1.21 The FWHM of the emission 
line was related to Mbh and Lb 0 i with a log-Gaussian 
distribution centered on the value given by Equation Q] 
with a scatter of ctfwhm (jShen et al.ll2008 ). To mimic 
the Mbh bias, we calculated the virial MBH.vir by using 
the distributions of LboI (converted to XL\ assuming a 
bolometric correction of 5) and dfwhm and Equation [0 
Note that the scatter of the virial Mbh to true Mbh 
is ^ (0.5ctt,) 2 + (2o~fwhm) 2 — 0.4 dex (IShen fc Kellvl 
120101) . We then apply the AGN luminosity and SFR 
cuts to the distribution of MBH,vir/M* to estimate the 
selection biases (i.e., a superposition of the Lauer et al. 
2007 and the M B h biases). 

We performed Monte Carlo simulations (based on the 
previous arguments) to estimate the selection biases in 
our sample. We started b y cons i dering; the galaxy stellar 
mass function of iMuzzin et al.l (120131) for star-forming 
galaxies: 

$(s) = (In 10)$*10 (s “ s * )(1+O!) exp(—10 s_s *). (10) 

a, s*, and 4)* (at each r e dshift bin) can be found from 
Table 1 of IMuzzin et al.l (I2013D (note that the absolute 
value of 4>* is not important for the bias estimation). For 
each M* that was randomly drawn from the distribution 
of Equation 1101 Mbh was randomly generated from the 
probability density function (PDF) o f Equatio n Qs] with 
ci = 1, c 2 = —2.85, and <r M = 0.3 (jHaring fc Rbdlfool) . 
In addition, for each M*, the SFR was randomly gen¬ 
erated from the PDF of Equation [5] with ap = 0.34 
(|Whitaker et al.l I2012D . The AGN luminosity and line 
width were calculated as described in the previous para¬ 
graph by assuming <tl = 0.4 and ctfwhm = 0.17 (such 
that the uncertainty of the virial Mbh is ~ 0.4 dex; we 
also tested other values of ctl and ctfwhm and found that 
the bias is not sensitive to <tl but depends more on the 
average value of the Eddington ratio distribution—see 
Section [4.1.21 for details). We repeated the simulation 
10 6 times to get the simulated data, i.e., our mock sam¬ 
ple. The PDF of these mock galaxies was estimated with 
kernel density estimation (green contour in Figure 0). 

We applied our sample selections to this mock sam¬ 
ple, starting with the limits on AGN luminosity and 
broad-line width. For AGNs in COSMOS, which form 
the bulk of our sample, t hese l imits are ?a b <22.5 and 
cfwhm > 10 3 kms -1 ([Trump et al.ll2009at ). The result¬ 
ing sub-sample is denoted as sample A. With these flux 
and wid t h limit s , our sample A is strongly affected by 
the lLauer et al.l d2007i) bias, as seen by the yellow con¬ 
tour in Figure 0 In the left panel (for 2 = 1), for 
example, we can clearly see a large offset between the 
Mbh — M* relation of our sample A and that of our 
mock sample. More quantitatively, (A(log Mbh/M*)) = 
(log(M BH /M*)) — c 2 = 0.2 for sample A. 

Next, a redshift-dependent cut-off in SFR was also 
applied to sample A to mimic the Herschel detection 
limit. The redshift-dependent SFR limit was estimated 
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from the flux limit of Herschel /PACS in COSMOS (since 
most of our sources are from the COSMOS field; the 
li mit is v L u (60 /mi) — 10 45 erg s _1 for z — 1, see, e.g., 
iRosario et al.ll2012h . The resulting sub-sample is denoted 
as sample B and is plotted as red points in Figure 0 
For the left panel, the bias ((A(log Mbh/M*)} = 0.11) 
is much smaller compared with sample A. We also ob¬ 
served similar features in the right panel (for z = 2): 
for sample A, (A(log Mbh/M*)} = 0.3; for sample B, 
(A(logM B H/M*)} = 0.22. 

We also compared the SFR — M* relation of sample B 
with the star-forming main sequence and found that 
sources in sample B show higher star formation activity. 
This result is simply due to selection effects: (1) there is 
intrinsic scatter in the “main-sequence” relation; (2) the 
galaxy stellar mass function is bottom heavy. That is, 
a star-forming galaxy with given SFR is more likely to 
be less massive, resulting in an apparent offset from the 
“main-sequence” relation. At z ~ 1, the median redshift 
for our sample, the (SFR/SFR ms ) of sample B is ~ 2: 
the same as the mean offset we obtained in Section [3~F1 
and Figure [6] This result suggests that, after accounting 
for selection biases, our sources are consistent with be¬ 
ing drawn from “main-sequence” galaxies. Therefore, it 
is valid to use the main-sequence relation in the simula¬ 
tions. 

For each redshift, we can perform similar simulations 
and therefore estimate the corresponding selection biases 
in sample B. Through this procedure, we estimated the 
selection biases of sample B as a function of redshift. 
We used this result to obtain the bias-corrected F1R04 
relation (i.e., a summation of the F1R04 relation and the 
bias). This bias-corrected relation will be included in 
Section 02] (see the red solid line in Figure (H). 

4.1.2. More Detailed Consideration of Biases 

The way we model the bias depends on the Mbh — M* 
relation, the main-sequence relation, their intrinsic un¬ 
certainties, and the Eddington-ratio distribution. As will 
be shown in Section [OJ our measured uncertainty in the 
Mbh — M* relation is ~ 0.3 dex, consistent with our as¬ 
sumed value. The uncertainty of the star-forming “main- 
sequence” relation is not prec isely known. If we ins tead 
use crp = 0.25 as reported in I Whitaker et al.1 (120121) for 
the “normal” star formation sequence, the bias of sam¬ 
ple B at e.g., z = 1 decreases to ~ 0.06 (the bias at 
z ~ 2.0 remains almost the same). We conservatively as¬ 
sume ctf = 0.34 when estimating the bias of our sample. 
The bias also depends on the connection between AGN 
activity and star formation. Let us assume, for exam¬ 
ple, the AGN luminosity is fully determined by galaxy 
properties (e.g., gas fueling) and is not directly linked 
with M B h- If so, there is n o bias for t he Mr h — A/* rela¬ 
tion (as already noted bv lLauer et alJ 12007 ). The reason 
is simple: in this case, our sample would be based on 
galaxy properties rather than M B h, just like the local 
sample. Some works suggest a li nk between AG N ac¬ 
tivity and galaxy properti es (e.g. iMu llaney et~akl 120121 
IRosario et all 120121 l2013allbl : iChen et al.l 1201 3L ~aiid this 
work). However, these AGN-galaxy relations may have 
large intrinsic scatter, perhaps making them inapplicable 
for modeling the bias in individual AGNs. 

We also have limited knowledge of the distri¬ 
bution of the Eddington ratio for BLAGNs (e.g., 


Kollmeier et al.l 2003 Shen et al. 

2008|; Trump et al. 

2009b| Hopkins et al.l 2009) 

Trump et al.l 201 it 


Kelly fc SlierJ 120131 b In our simulations, we also 


tried other Eddington-ratio distributions, e.g., a log- 
uniform distr ibutio n (i. e., close t o the “fiducial” model 
proposed by iHickox et al.l 12014] ) or a log-normal dis¬ 
tribution with a scatter of 0.8 dex. We additionally 
required these distributions to be truncated at both 0.01 
and 1.0. These two cut-offs are m otivated by bot h the 
theory of accretion dis ks (see, e.g, iNaravan fc Yi 1995; 
Yuan fc Naravanll2014fl and observations (e.g.. IH 0 II 2 OOS : 
Trump et al.l I2009bl 1201 If) . We found that all these 
distributions give similar results regarding selection 
biases as long as the mean value of these distributions 
is the same (i.e., 0.1). If we instead assume a smaller 
mean value for the Eddington ratio, the selection biases 
increase (and vice versa if we assume a larger mean value 
for Eddington ratio). This is simply due to the fact that 
the selection limit of Mbh increases with decreasing 
Eddington ratio and the slope of the Mbh function 
becomes steeper with increasing Mbh- To illustrate 
its effect, we adopted a power-law distri butio n of the 
Eddi ngton ratio with a slope of —1.0 ()Bongiorno et al.1 
120121 1 and a cut-off of > 0.01 and < 1.0. This distri¬ 
bution has a mean Eddington ratio of ~ 0.02. The 
resulting selection bias on the M B h — M* relation 
is 0.22 (0.32) at z = 1 (z = 2). Future constraints 
on the Eddington ratio distribution are crucial for a 
more reliable estimation of the biases. Regardless, we 
emphasize that our conclusions in Section 14.21 would not 
be significantly changed if the selection bias is essentially 
zero or as large as ~ 0.3. 

The detectability of a BLAGN actually depends on 
the luminosity contrast between the AGN and the host 
galaxy, since this determines whether we can detect 
the broad H/3 and/or Mgli emission lines (and there¬ 
fore measure Mbh)- This effec t has not be en discussed 
in previous simulations of the lLauer et al.l (l2007f) bias. 
To check the luminosity contrast in ou r sim ul ated sam¬ 
ple s, we used t he 3D-HST catalogs (iBrammer et al.lf2012b 
iSkelton et al.1 120141) P°1 The 3D-HST catalogs have red¬ 
shift, stellar mass, age, SFR, and rest-frame colors for 
inactive galaxies and can be used to assign galaxy prop¬ 
erties to our simulated samples. For each redshift bin, we 
started by dividing our mock sample into galaxy stellar 
mass bins. As a second step, for each redshift and galaxy 
stellar mass bin, we selected main-sequence galaxies (de¬ 
fined by Equation [7] and its 3cr uncertainty) from the 3D- 
HST catalogs whose redshifts and galaxy stellar masses 
lie within the bin. Then we used the selected sources to 
determine the distributions of the galaxy luminosity at 
(rest-frame) 2800 A and the SDSS-g band (whose effec¬ 
tive wavelength is 4718.8 A, close to the wavelength of 
H/3) for each bin. Finally, we assigned the (rest-frame) 
galaxy luminosity to our mock samples according to the 
distributions. We again applied the limits in AGN lu¬ 
minosity, FWHM, and SFR. For the resulting samples, 
we calculated the distribution of the luminosity ratio of 
the AGN to the host galaxy at 2800 A (for sources at 
z > 0.88, whose M B h are measured using Mgli) or the 
SDSS-g band (for sources at 2 < 0.88 whose Mbh are 

20 http://3dhst.research.yale.edu/Home.html 
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measured using H/3). We found that for the mock sample 
A, there is a fraction of sources with Lagn/-^G ai < 0.1 
at these wavelengths. However, no source in the mock 
sample B (i.e., the one which suffers the same selection 
biases as our data) has a luminosity ratio < 0.1. This 
suggests that our assumption of the AGN contribution 
to the total SED in Section EO is appropriate and does 
not introduce additional selection biases to our results. 

In order to measure the galaxy total stellar mass via 
two-component SED fitting (see Section lT2l) . the galaxy 
must be detectable and not outshone by the AGN in 
the K band. This requirement will introduce an upper 
limit to Mbh/M*. We calculated this limit for the galaxy 
emission to be at least 10% of the AGN emission in the K 
band. Assumin g an AGN bolometric correction of ~ 10 
at the K band dRichards et, al.ll2006D and an Eddington 
ratio of 0.1 (see Fi gure Band e g.. iKollmeier et alJl200fit 
iTrumn et al.ll201 it ILusso et al.ll2012h . the A-band AGN 
luminosity is 

Aagn = 1.3 x 10 36 Mbh. (11) 


Meanwhile, the i^-band galaxy luminosity can be written 
as 


K _ ^Gal T 
^Gal — -k©> 


( 12 ) 


where 7 is the mass-to-light ratio. We ad opted log(y) ~ 
—0.6 for our z ~ 1 star-forming galaxies (lArnouts et al.1 
120071 with a correction for the IMF of Chabrier 2003). 
Note that there is a weak dependence of 7 on redshift, 
although this is not significant over the redshift range of 
our sources. With these parameters the requirement of 
L Gal > 0-1-^AGN indicates that Mbh/M* cannot exceed 
0.16. This Mbh/M* limit is more than two orders of 
magnitude larger than the mean value of the Mbh ~ M* 
relation (which suggests (Mbh/M*) ~ 0.001). Therefore, 
we expect this limit should not introduce significant bi¬ 
ases to our data. We also directly test the effect of the 
L Gal > 0.1L^ GN limit by assigning the rest-frame K- 
band luminosity to our Monte Carlo simulations of the 
Mbh — M* relation, again finding that the galaxy detec¬ 
tion limit does not introduce significant bias. 


4.2. Black Hole Mass-Galaxy Total Stellar Mass 
Relation of Star-forming Galaxies 

In this subsection, we explore the Mbh — M* rela¬ 
tion of BLAGNs in star-forming galaxies. The upper-left 
panel of Figure |5] plots the distribution of our sources 
in the SMBH mass-galaxy total stellar mass plane. For 
comparison, we also include the sample of local inac¬ 
tive galaxies from HR04 and the local SMBH mass-bulge 
mass relation of HRO^EH and its la uncertainty. As 


21 iKormendv fc Hoi d2013l f point out that the HR04 SMBH-to- 
bulge mass ratio is probably underestimated by a factor of ~ 2 
(see Section 6.6 in their paper). However, our single epoch virial 
SMBH mas s estimators are ca librated by the local Mbh — n rela¬ 
tion of ITremaine et all (2002) which is consistent with the HR 04 
SMBH mass-bulge mass relation. That is, the geometry factor 
which converts the virial product into a SMBH mass is derived 
from a Mbh — a relation that is consistent with the HR04 relation 
IQnken et al.ll2O04 ). Therefore we retain the older HR04 relation 
for consistency wi th the single epoch viria l SMBH mass estimates. 
Updating to the IKormendv K: Hoi ( 2013 ) relation would require 
similar updating single epoch virial SMBH estimates by (roughly) 
the same factor, resulting in no difference to our conclusions re¬ 
garding evolution. 


seen from the upper-left panel of Figure [8j our sample 
is consistent with the HR04 relation. In fact, only 13/69 
sources lie outside the lcr uncertainty of the HR04 rela¬ 
tion when considering error bars. A Mann-Whitney U 
test of Mbh/M* on our sample and the local inactive 
galaxies also indicates that the two samples are statisti¬ 
cally consistent (with a null probability p = 0.2). The 
same conclusion holds for our low-redshift sub-sample, 
although considering the high-redshift subsample alone 
results in a small (but statistically significant) difference 
from the HR04 relation. However, this is likely due to 
selection bias since the bias increases with redshift. In¬ 
deed, if we take the selection bias (see Section 14.11) into 
consideration, the Mann-Whitney U test indicates that 
there is no significant difference between the high-redshift 
sub-sample and the local sample of HR04 (with p = 0.5). 
We also binned our sources by M*. Four bins were cre¬ 
ated: M* < 10 10 - 5 M q 10 10 ' 5 Mq < M* < 10 11 M 0 , 
10 11 Mq < M* < to 11 - 3 Mq, and M* > 10 115 M 0 . We 
calculated the median values of Mbh and M* and their 
lcr uncertainties for each bin. The upper-right panel of 
Figure [8] plots the results. As a reference, we again in¬ 
clude the sample of local inactive galaxies from HR04 
and the local SMBH mass-bulge mass relation of HR04 
and its lcr uncertainty. As seen from this panel, the 
binned data are consistent with the local SMBH mass- 
bulge mass relation of HR04. 

We can also quantify possible evolution by testing the 
correlation between Mbh/M* and redshift. The lower 
panel of Figure [8] plots this mass ratio and its uncer¬ 
tainty as a function of redshift. As a reference, we also 
plot the local mass ratio (black dashed line) and its lcr 
uncertainty. For each redshift, we estimated the cor¬ 
responding selection bias in Section 14.11 The red solid 
line in the lower panel of Figure [8] shows the expected 
HR04 relation after considering the selection bias. We 
also calculated the mean value of Mbh/M* and its 2cr 
uncertainty for sources with z < 1, 1 < z < 1.5, and 
z > 1.5; red squares show the result. Comparing the red 
squares with the red solid line, we find our sample shares 
the same mean value of Mbh/M* with the bias-corrected 
local one. 

We also used the simulation-ext rapolation (SIMEX) 
technique dCook fe StefanskH I1994D to estimate the in¬ 
trinsic scatter (subtracting the measurement errors) of 
the distribution of Mbh/M*. The SIMEX estimation 
gives a scatter of 0.36 dex for our sample. This indicates 
that the intrinsic scatter of Mbh/M* for our sample is 
not significantly larger than the local value. Our ob¬ 
served scatter in Mbh/M* is somewhat affected by the 
biases, so we do not put too much emphasis on the sim¬ 
ilarity in scatter. Still, the similar scatter suggests that 
our estimation of selection bias is reliable, and that our 
sample is consistent with the HR04 relation (also see Sec- 
tion l6.1D . This non- e volution result agrees we l l with , e.g., 
iJahnke et al.l (|2009D : iSchramm Sz Silvermanl (120131). but 
it con trasts with some other works, e.g.. iMerloni et al.l 
(2010). The latter claimed an evolution of Mbh/M* 
at the > 5 a significance level, although this is proba¬ 
bly only due to the select ion biases (jMerloni et alJl20lfj 
iSchulze fe Wisotzkill2014D . 
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5. TIME EVOLUTION OF THE SMBH 
MASS-GALAXY TOTAL STELLAR MASS 
RELATION 

In the previous section, we find that the Mbh — A/* 
relation does not evolve with z up to z ~ 2. In addi¬ 
tion, the X-ray and FIR data allow us to calculate the 
current mass growth rate of both AGNs and their host 
galaxies. We use these mass growth rates to assess if our 
galaxies and AGNs will remain on or move away from 
the established Mbh — M* relation. 

We start by exploring the distribution of the specific 
SMBH mass growth rate (defined as s M = M/Mbh) 
and the specific galaxy total stellar mass growth rate 
(defined as sSFR = SFR/M*). Figure [9] plots the spe¬ 
cific mass growth rate distributions for the high-redshift 
and low-redshift sub-samples. There are two points to 
emphasize: (1) in both sub-samples, the instantaneous 
specific SMBH growth rate is not smaller than that of 
the galaxy (especially for the low-redshift sub-sample, 
the instantaneous specific growth rate of the SMBH is 
much larger than that of the host galaxy); (2) the spe¬ 
cific SMBH growth rate does not evolve much (the mean 
offset is ~ 0.2 dex) while the specific galaxy growth rate 
decreases significantly (the mean offset is ~ 0.8 dex) from 
high to low redshift. That is, our results indicate that 
the apparent sM/sSFR seems to increase from high to 
low redshift. This turns out to be largely due to selec¬ 
tion effects; after modeling the biases in Section [Gl we 
find that the data are consistent with a uniform AGN 
duty cycle at both high and low redshift. 

In addition, for both sub-samples, the specific black 
hole mass growth rate distributions show negligible tails 
at the growth rate < 0.2(10~°') Gyr -1 ~ O.OlMEdd (see 
Figure where M^dd = ^Edd/OTc 2 . T hese neg ligi¬ 
ble ta ils agree with the c onclus io n of, e.g., iCollin et al.1 
120061): i Kollmei er et al.1 (|200fiH : I Trump et all (120 1 ID : 
ILusso et all (I2012h that BLAGNs have a minimum Ed¬ 
dington ratio of ~ 0.01. 

We use these specific growth rates to i nvesti gate the 
“flow patterns” (following iMerloni et al.1120101) of our 
sample in the Mbh — M* plane. Figures [TU] and HH 
plot the results. The direction of the vector is defined 
by sM/sSFR while the color map indicates the absolute 
specific growth rate (i.e., a summation of the specific 
SMBH mass growth rate and the specific galaxy total 
stellar mass growth rate in quadrature) along the vec¬ 
tor (which spans more than two orders of magnitude). 
We again include the local HR04 relation and its lcr un¬ 
certainty as the dashed line and the shaded region. For 
simplicity, these figures do not show the uncertainties of 
sM/sSFR on the plot (the uncertainties can instead be 
seen in Figure IT2l) . In general, AGNs and galaxies which 
are outliers in Mbh/M* tend to have evolutionary vectors 
whose directions are anti-correlated with their positions: 
galaxies with over-massive SMBHs tend to have lower 
sM/sSFR, and galaxies with under-massive SMBHs tend 
to have higher sM/sSFR. This anti-correlation behavior 
can be seen not only for the whole sample (Figure fTTIl) 
but also for each sub-sample (the left and right panels 
of Figure fTTl) . As will be shown in Section 1B~21 this re¬ 
sult is not solely due to selection effects. In Sections [6.21 
and l6.3l we demonstrate the anti-correlation between the 


(instantaneous) growth rate and mass ratios, if extended 
over some time with an AGN duty cycle of 0.1, implies 
that AGN activity and star formation could maintain the 
Mbh — M* relation. 

6 . DISCUSSION 

6.1. Building the Local Mbh — Mb u i Relation 

In Section 14.21 we showed that our sample exhibits no 
apparent redshift evolution in the Mbh — M» r elation. 
There fore, our results along with , e.g., IJahnke et al.1 
(120091) and ISchramm fe Silverman! (120131 1 suggest that 
the Mbh — M* relation does not evolve strongly with 
redshift. As stated before, the local inactive galaxies of 
HR04 are bulge-dominated (Mb u i — M*). Meanwhile, 
our sources (at least a significant fraction) are likely to 
have substantial disk components. First, our sources 
show strong star formation activity which indicates non- 
negligible disk components (for example, Lang et al. 
2014 suggest that star-forming galaxies are likely to have 
smaller bulge-to-total ratios). In addition, for the X-ray 
luminosity and redshift ranges of our sources, a substan¬ 
tial fraction o f AGNs reside in d i sk galaxies (i.e., Mr„i <C 
M«; see, e.g.. IGabor et al.ll2009t IGeorgakakis et al.ll2009l 

iCisternas et al.1 120111 : iKocevski et all 120 1 2fu Further¬ 

more, we have als o cros s-matched our sample with 
that of iGabor et aD (120091) (for COSMOS sources) and 
ISchawinski et al.1 (120111 ) (for CDF-S sources) and found 
that five out of the fourteen matched sources have Sersic 
index much smaller than 2.5. Therefore, the Mbh — Mb u i 
relation is likely to evol ve (at l east mildly) with redshift. 
As suggested bv IJahnke et all (l2009h . a redistribution of 
stellar mass from disk to bulge is required to build the 
local Mbh — Mb u i relation. 

Our results suggest that in the next ~ 10 Gyr a typical 
galaxy in our sample would (1) increase its galaxy total 
stellar and SMBH masses accordingly; (2) transfer much 
of its stellar mass from its disk to its bulge. The former 
indicates that gas fueling controls both AGN activity and 
star formation (see Section ITT~2l for more details). For the 
redistribution of stellar mass from disk to bulge, possi¬ 
ble mechanisms are, e.g., violent disk instabilities (which 
are driven by co ld gas, see e.g ., Dckcl et al. 2009|) and/or 
galaxy mergers (iHopkins et al.1 20101) . The former mech¬ 
anism is proposed since our sources are likely to have sub - 
stantial cold gas (jRosario et alJl2013at iVito et al.ll2014[) . 
although observational evidence of violent disk insta¬ 
bilities in AGN ho sts is mixed (jBournaud et al.1 120121 
iTrump et al.l [20141. The latter mechanism is also at¬ 
tractive as it can not only redistribute stellar mass 
from disk into bulge but also help reduce the i ntrin¬ 
sic sc atter in the Mbh — M* relation (.lahnk e fc Maccibl 
2011). Note that our AGNs are more likely fueled by 
cold gas rat her than triggered by galaxy-galaxy merg¬ 
ers (see, e.g.. IGabor et al. 120091 IGeorgakakis et al .1120091 
ICisternas et al.ll201lHKocevski et al.ll2012l and our next 

Section). Therefore, if galaxy-galaxy mergers do play a 
role, they are likely to be dry mergers rather than gas- 
rich wet mergers. 

Currently, our estimation (see Section 14.21) indicates 
that the intrinsic scatter of the Mbh — M* relation is not 
significantly larger than that of the HR04 relation. How¬ 
ever, this estimation is also subject to selection biases. 
Using our simulated samples in Section 14.11 we found 
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that a range of intrinsic scatter (e.g., 0.36— 0.5 dex) can 
reproduce the observed scatter. On the other hand, our 
selection corrections in Section [472] are made by assuming 
the Mbh — M* relation has the same intrinsic scatter as 
the HR04 relation and our data are consistent with the 
bias-corrected HR04 relation. This indicates that the in¬ 
trinsic scatter of our sample should not be significantly 
larger than that of the HR04 relation (i.e., the intrin¬ 
sic scatter is at most 0.5 dex). An average of one major 
merger for each of our galaxies is enough to reduce the in¬ 
trinsic scatter by a factor of y/2 (IJahnke fe Macci6ll2011fl 
if the intrinsic scatter is indeed 0.5 dex rather than ~ 0.3 
dex. 

6.2. Growth of SMBHs and Host Galaxies 

In Section [5j we demonstrated that AGNs which are 
offset from the Mbh — M* relation tend to have (instan¬ 
taneous) evolutionary vectors which are anti-correlated 
with their mass offsets. That is, galaxies with over- 
massive (under-massive) SMBHs tend to have lower 
(higher) sM/sSFR. We quantify the anti-correlation be¬ 
tween Mbh/M* and sM/sSFR by performing a linear fit 
(via chi-squared minimization, taking the uncertainties 
into consideration). Results are plotted in the upper pan¬ 
els of Figure [12] The red vertical dashed lines represent 
the local Mbh/M*. The red horizontal dashed lines cor¬ 
respond to s M /sSFR = 1. There are more sources above 
the horizontal lines. If the growth rate rates persist over 
some length of time, then this asymmetric distribution 
indicates the AGN duty cycle is less than unity (i.e., the 
lifetime of an active SMBH is smaller than the lifetime of 
star formation); otherwise the SMBH will become over- 
massive as we evolve our sources. For example, if we 
assume our sources keep their current M and SFR for 
some timescales £agn and fsF, we expect BLAGN hosts 

to evolve as A (log M B h/M*) = log( 1+^^ )■ To main¬ 
tain the Mbh — M* relation, we expect that, on average, 
£agn < tsF- In Section 16.31 we derive the appropriate 
AGN duty cycle (^agn/Isf) such that, under appropriate 
assumptions about the AGN activity and star formation 
histories, the Mbh — M* relation is maintained. The red 
solid lines are for the best fits of the anti-correlation: 
for the high-redshift sub-sample, the best fit is y = 
—2.7(±0.3) — 1.14(±0.12):e; for the low-redshift sub¬ 
sample, the best fit is y = —1.7(±0.4) — 0.86(±0.12):c. In 
each sub-sample, the slopes are consistent with —1 and 
differ from zero at a high significance level (for the high- 
redshift sub-sample, the null probability is p = 4 x 10 -12 ; 
for the low-redshift one, p = 2 x 10~ 7 ). 

The anti-correlation between the specific growth rate 
ratio (sM/sSFR) and the mass ratio, with slope of —1, 
implies that the absolute growth rate ratio (M/SFR) is 
not correlated with the mass ratio and is a constant (with 
small scatter). We demonstrate this in the lower pan¬ 
els of Figure 1121 A linear fit (via chi-squared minimiza¬ 
tion) between M/SFR and Mbh/M* results in a slope 
of 0.27 ± 0.3, i.e., statistically consistent with zero. The 
relatively narro w distribution o f the M /SFR ratio is in 
agreement with iMullanev et al.l (12012D . after accounting 
for the AGN duty cycle (see Section [fOll . 

To test whether this anti-correlation and the narrow 


distribution of M/SFR ratio are simply due to selection 
effects we again turn to the Monte Carlo simulations in 
Section [4] These simulations created mock samples with 
Mbh, M*, M, and SFR following the same selection cri¬ 
teria as our data. We tested three Eddington-ratio dis¬ 
tributions which are log-normal with scatter (ox) of 0 . 8 , 
0.6 and 0.4, respectively. Meanwhile, we use the same 
SFR — M* and Mbh — M* relations and scatter as in 
Section 4, and ctfwhm is set such that the uncertainty of 
the virial A7 bh is — 0.4 dex. We created 69 mock sam¬ 
ples, one at each redshift of the 69 observed BLAGNs, 
with 4 x 10 5 sources in each mock sample. 

The initial mock sample (i.e., the one without any se¬ 
lection cut-offs) at each redshift does not show an anti¬ 
correlation between the ratio of the masses and their 
growth rates. Indeed, Spearman’s correlation test in¬ 
dicates a coefficient of p ~ —0.002 for ctl = 0.4 or 
weaker depending on ox). We then applied the sam e 
(redshift-dependent) cut-offs presented in Section |4~T1 to 
these mock samples. Finally, we randomly selected 69 
sources from the 69 mock samples (i.e., selected one 
source per mock sample) and calculated the correspond¬ 
ing p (for our BLAGNs sample, p = —0.78) and the 
scatter of the distribution of M/SFR (for our BLAGNs 
sample, the scatter is 0.4). Note that, in this step, we 
added 0.3 dex uncertainty to M and 0.17 dex uncertainty 
to SFR to mimic approximately the measurement error 
in each variable (see Section 1531) . We repeated the se¬ 
lection 10 6 times and constructed the distribution of p. 
We found that, of all the Eddington-ratio distributions 
tested, the ox = 0.4 model comes closest to reproduc¬ 
ing the observed anti-correlation between s M /sSFR and 
Mbh —M*, albeit with a probability of only ~ 0.5%. The 
simulations also indicate that the probability to repro¬ 
duce the observed scatter of the distribution of M/SFR 
is only ~ 0.05%. Note that in the ox = 0.4 model, 
the scatter between AGN luminosity and SFR is only 
~ 0.6 dex. Broader Eddington-ratio distributions (e.g., 
uniform or shallow power-law) decrease the likelihood of 
selection effects reproducing a constant M/SFR ratio, as 
do larger values of oy, or of. Moreover, we also performed 
the same simulations to the high- and low-redshift sub¬ 
samples and found the likelihood of selection effects also 
decrease. 

From the simulations, we can conclude that the ob¬ 
served anti-correlation between specific growth ratio and 
mass ratio and/or the narrow distribution of M/SFR 
are unlikely to be solely due to selection effects but in¬ 
stead support the idea that there is a (instantaneous) 
connection between star formation and SMBH accretion 
activity (see the lower panels of Figure m- However, 
this conclusion rests upon the (rather uncertain) assumed 
Eddington-ratio distribution, as well as the similarly un¬ 
certain intrinsic scatter in Mbh and Mbh — M*. Better 
estimates of these quantities and their scatter are nec¬ 
essary to fully confirm that our results are not purely a 
selection effect. 

AGN feedback may explain the anti-correlation be¬ 
tween sM/sSFR and Mbh/M, and the constant ratio of 
M /SFR. As has been long suggested, AGN feedback may 
regulate star formation and therefore maintain a corre¬ 
lation between AGN accretion and galaxy star forma- 
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tion ( e.g.-lSilk fe Reedll998t lKiridl20()3l : iDi Matteo et al.1 
120051: iHopkins et al.l l2006f) . However, recent simula¬ 


tions with a more realistic ISM structure suggest that 
AGN feedback can inject substantial energy into the 
ISM without h aving significant effect on cold gas in 


the galaxy (e.g., 

3ourne et al.112014 Gabor & Bournaud 

2014: Roos et al. 

20TD- If so, our results are difficult to 


explain by AGN feedback. 

Our results can also be understood in a scenario 
where AGN activity and star formation are both gov¬ 
erned by cold gas supply, without a need for AGN 
feedback. This scenario could lead to a connec¬ 
tion between the AGN luminosi ty and SFR or even 
an AGN “main sequ en ce” (e.g.. IMullanev et al.l l2012t 
iRosario et al.1 l2013allbl : lYito et al.l 120141) . For a, host 
galaxy with an under-massive SMBH, we would expect 
the specific SMBH mass growth rate to be larger than 
that of the host galaxy since the mean ratio of the SFR to 
AGN lumin osity remains roughly constant with changing 
M* (IMullanev et al.ll2012f) . That is, an under-massive 
SMBH has a preferentially steeper evolutionary vector 
than an over-massive SMBH on the Mbh — M* plane. 

6.3. The AGN Duty Cycle among Star-Forming 
Galaxies 

We are also interested in the AGN duty cycle and 
star formation timescale, tsf, that would maintain the 
Mbh — M* relation as our sources evolve forward in 
time. We started with a Monte Carlo simulation which 
re-sampled Mbh, M*, M, and SFR based on the ob¬ 
served values and errors for every source in our sample. 
For each realization, we evolved every source in time 
and calculated their new positions in the Mbh — M* 
plane by integrating SFR(t) and M(f) over a grid of 
AGN duty cycles and tsf- For the high- (low-) red- 
shift sub-sample, all sources were evolved from their 
current redshifts to z = 1 (z = 0.2). SFR(f) was as¬ 
sumed to be SFR(f) = SFR(i 0 ) exp(— (t— to)/ T SF), where 
SFR(f 0 ) was the re-sampled SFR. As for M(t), we as¬ 
sumed that, on average (e.g., over ~ 100 Myr), M{t) = 
M(to) exp(— (t — Io)/tagn)- That is, following similar 
derivations in the previous section, the new positions of 
our sources are (M B H,new, M* jnew ), where M B H,new = 
M B h ,init H - f M(t)dt ~ MBH.init + ^agnM^o) and 
M*, new = M* 

,init / SFR (t)dt ~ M* ,init + Ts F SFR(f 0 ) 
(as long as the evolution timescale is much larger than 
tsf)- MBH.init and M * t i n it are the observed SMBH mass 
and galaxy stellar mass, respectively. We define the AGN 
duty cycle as tagn/tsf- Assuming our AGN hosts are 
drawn fr om the star-forming main sequence —as sug¬ 
gested bv lRosario et al.1 (j2013bl ) — this definition of AGN 
duty cycle essentially measures the fraction of BLAGNs 
among star-forming galaxies (we will discuss this defini¬ 
tion later). We then determined the fraction (/ on ) of our 
sources whose new positions are consistent with the bi¬ 
ased HR04 mass ratio (see Section ITU and its 1 <t uncer¬ 
tainty for each combination of AGN duty cycle and tsf- 
Finally, we found the AGN duty cycle and tsf that result 
in the highest f on . After 10 4 re-sampled realizations, we 
obtained the distribution of the AGN duty cycle and tsf 
that would best maintain the Mbh — M* relation. 

The duty cycle used here is only equivalent to 


sSFR/sM or SFR/M if s M and sSFR are both very 
small (i.e., if sMtagn and sSFRtsf both are <C 1 , then, 
A(logM B H/M*) « 1ii(10)(sMtagn — sSFRtsf))- This is 
true for our low-redsliift sample. If, however, sMtagn 
and/or sSFRtsf are ~ 1 (e.g., our high-redshift sample), 
the AGN duty cycle is a complicated function of Mbh, 
M*, M, and SFR. Actually, for our assumed SMBH ac¬ 
cretion and star formation histories, our definition of the 
duty cycle measures the ratio of the lifetime of BLAGNs 
to that of active star formation. 

The two-dimensional PDF of the AGN duty cycle and 
tsf is plotted in Figure [13] for the high-redshift (left 
panel) and low-redshift (right panel) sub-samples. The 
blue contour in each panel represents the “la” scatter of 
the most likely AGN duty cycle and tsf- That is, given 
the measurement errors of Mbh, M*, M, and SFR, there 
is a 68.3% probability that AGN duty cycles and tsf 
which best maintain the Mbh — M* relation are within 
the contour. To estimate the one-dimensional distri¬ 
bution of the AGN duty cycle (regardless of tsf), we 
integrated the two-dimensional PDF along tsf- From 
this distributio n, the preferred AGN duty cycle is ~ 0.1 
(agreeing with iSilverman et al.l 120091 ) in both the high- 
and low-redshift sub-samples (O.lloios at 2 > 1 , and 
O.ll 0 : 0 9 at 2 < 1 , with errors defined by the 15.87 and 
84.13 percentiles). In other words, our data favor a 
non-evolving (by no more than a factor of ~ 4) AGN 
duty cycle of ~ 10 % for star-forming galaxies at both 
z ~ 1.5 and z ~ 0.5, ruling out very high (~ 67.6%) 
duty cycles at the 95% confidence level. The large con¬ 
fidence intervals for the AGN duty cycle are partly due 
to the large uncertainties in the estimated masses and 
growth rates. There may also be a large intrinsic scatter 
in AGN duty cycle between different galaxies. We note 
that these duty cycles are appropriate only for rapidly ac¬ 
creting SMBHs (i.e., BLAGNs or other AGNs with sim¬ 
ilar Eddington ratios) in star-forming galaxies. Lower- 
luminosity AGNs are likely to have higher duty cycles, 
while quie scent galaxies (known t o have lower AGN f rac- 
tions, e.g. ■ IRosario et al.ll2013al lbl: lTrump et al.ll2013f) are 
lik ely to have lower duty cycles (for example, as found 
bv lBongiorno et aLll2012i) . 

7. SUMMARY AND FUTURE WORK 
7.1. Summary 

In this paper, we studied 69 BLAGNs selected from the 
COSMOS and CDF-S fields based on X-ray observations, 
FIR observations, high-quality spectroscopy, and multi¬ 
band photometry. With these data, we simultaneously 
determined Mbh, M*, M, and SFR and investigated the 
evolution of the Mbh ~ M* relation. Our main conclu¬ 
sions are the following: 

(1) Our sample suggests that, up to z ~ 2.0, there is 
no evolution (no more than ~ 0.2 dex) of the Mbh — M* 
relation. The Mbh — Mb u i relation, however, is likely 
to evolve, simply because our galaxies are expected to 
be more disk-dominated than the local sample of HR04. 
See Sections 14.21 & 16.11 

(2) Our data indicate an anti-correlation between 
Mbh /M* and the ratio of the two specific mass growth 
rates. Such an anti-correlation suggests that the on-going 
AGN activity and host-galaxy star formation are phys- 
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ically connected. Under appropriate assumptions about 
the AGN accretion and galaxy star formation histories, 
this anti-correlation can maintain the Mbh/M* relation. 
See Sections ED&EO] 

(3) We also investigate the possible values of AGN duty 
cycle which best maintain the non-evolving Mbh ~ A/* 
relation. Our data favor a non-evolving (i.e., a factor of 
< 4) AGN duty cycle of about 0.1 for rapidly accreting 
SMBHs in star-forming galaxies (O.l^os at z > 1, and 

0 -llo’o9 at 2 < 1). See Section RTTTT1 

These results fit into a picture where the same gas 
reservoir fuels both AGN activity and galactic star for¬ 
mation. 

7.2. Future Work 

Our results can be advanced in several regards. First, 
our sample size is significantly limited by the Herschel 
sensitivity. Only ~ 1/3 of BLAGNs are detected by Her¬ 
schel. In this work we attempt to correct for this bias, 
but it would be better to simply have a less-biased sam¬ 
ple. Future deeper rest-frame FIR observations (e.g., 
with ALMA, SCUBA-2, SPICA) can enlarge the sam¬ 
ple size and are crucial for reducing the uncertainty and 
improving our understanding of the connection between 
AGN activity and star formation. As we have pointed 
out in Section [4Tl the biases we mentioned (e.g., the bias 
of Mbh alone and the bias of Mbh/M*) depend on this 
connection. Therefore, deeper FIR surveys can also help 
us better address the selection biases in this work. This 
improvement, in turn, providing benefits to both our un¬ 
derstanding of the co-evolution of SMBHs and their host 
galaxies (e.g., the formation of the local Mbh — Mb u i 
relation) and SMBH mass demographics. 

Also, the large intrinsic scatter of Mbh prevents us 
from obtaining higher confidence level conclusions. Fu¬ 
ture improvements of the single-epoch virial estimators 
will be very helpful (for example, with la rge multi-object 
reverberation mapping campaigns iShen et al.lf2014h . 

Finally, more detailed observations and studies of phys¬ 
ical properties (other than M*) of the host galaxies 
are also very important. These properties (e.g., mor¬ 
phology and gas content) can help place additional dis¬ 
tinct constraints on the co-evolution path of SMBHs and 
their host galaxies (e.g, the bulge-to-total ratio and the 
AGN/star formation triggering mechanism). 
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Fig. 1.— Examples of multi-component fits to the broad H/3 (the upper-left panel) or Mgn emission lines. Dotted lines represent minor 
features (e.g., the [O ill] AA4959, 5007 lines and the narrow H/3 component which are removed in the broad-line fit). 
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Observed-Frame Wavelength [A] 


Fig. 2. — Examples of two-component (AGN+galaxy) fits to SEDs (from the NUV to NIR). Red dotted (green dashed) lines represent 
the AGN (galaxy) component. Black solid lines show the summation of the galaxy and AGN SEDs. Note that each data point has a 
photometric error which is typically too small to be visible on this scale. 
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Fig. 3. — The bolometric luminosity of AGNs (calculated by combining the X-ray observations and SED fitting) vs. redshift. Blue 
triangles represent AGNs that are detected in the observed-frame hard X-ray band (2 — 10 keV); red points are sources that are detected in 
the observed-frame soft X-ray band (0.5 — 2 keV) but not in the observed-frame hard X-ray band. The green solid line repre sent the limi t 
introd uced by ^ab =22.5 which corresponds to the ~ 90% complete ness of the optical sp ectroscopy for COSMOS sources (1 Trump et al.l 
2009a). The bolometric correction of the i band is assumed to be 12 ( Richards et al. 2006) and the K correction is made by assuming the 
optical-to-UV SED is a power-law with a u = —0.44 JVanden Berk et al.H200U) . The three green dashed lines represent the limits introduced 
by the X-ray sensitivity of t he soft band of XMM -Newton /COSMOS and Chandra /COSMOS and the hard band of CDF-S, respectively. 
The bolometric correction of I Hopkins et all ^20071) is assumed to calculate LboI from the X-ray sensitivity. 
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Fig. 4.— The bolometric luminosity of AGNs vs. Mbh- The pink dashed, green solid and red dotted lines correspond to the Eddington 
ratios of 1.0, 0.1 and 0.01, respectively. For each object, the intrinsic scatter in the Lb 0 i and Mbh estimators (0.35 an d 0.4 dex, respectively) 
dominates the error. The Eddington ratios of our source spans from 0.01 to 1.0, typical range of BLAGNs (e.g., Kollmeicr et al. 2006; 
I Trump et a l. 2009b, 2011; lLusso et aLll2012l) . 
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Fig. 5. — Illustration of the reliability of estimating FIR luminosity using different Herschel bands. Left: Comparison of FIR luminosity 
estimated from two bands with that from all three bands. Right: Comparison of FIR luminosity estimated from one band with that from 
all three bands. The solid line and shaded area represent the ratio of unity and 0.17 dex uncertainties. The two-band FIR estimator agrees 
well with the three-band estimator. Estimating SFR from a single band is also consistent with the three-band estimator. Note that there 
is a source in the right panel shows an extremely large uncertainty which is caused by the fact that the uncertainty of the flux in 250 ^m 
for this source is large (if we take also the 100 ^m and 160 ^m bands into consideration, the uncertainty of this source is small). 
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Fig. 6.— A comparison of SFR we measured with the expected SFR of star-forming main sequence galaxies with the same stellar mass. 
Our sources have a star formation activity which is marginally higher than that of the “main- sequ ence” (the mean of SFR/SFR ms is 2). 
This offset is due to an Eddington bias which is driven by the Herschel flux limit (see Section Id.ll for more details). The red dashed line 
represents the bias corrected SFR/SFRms- 
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Fig. 7.— Monte Carlo simulation of biases in our data. Left (right) panel is for z = 1.0 (z = 2.0). The green and yellow contours 
correspond to the PDFs of our mock sample and sample A (which results from the mock sample with a AGN luminosity cut off). For sample 
B, which results from sample A with a SFR cut off (i.e., sample B and our data share the same selection bias.), the actual distrib ution of 
simulated sources in the Mbh — M* plane is plotted. The solid green lines correspond to log Mbh = logM* — 2.85. See Section |4.1| for 
more details. 
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Fig. 8.— Upper-left panel: The distribution of our BLAGNs in the Mbh — M* plane. Upper-right panel: Sources are binned by M*; 
four orange points represent the median values of Mbh and M* and their la uncertainties (estimated via bootstrapping) for sources within 
four bins: M* < 10 10,5 M©, 10 10,5 M© < M* < 10 11 M©, 10 11 M© < M* < 10 11,5 M©, and M* > 10 11,5 M©. The triangles in both 
panels represent sources from HR04. Lower panel: Logarithm of the ratio of the SMBH mass to the galaxy total stellar mass vs. redshift. 
Three red squares represent the mean value of the mass ratio and its 2 a uncertainty for sources with z < 1.0, 1.0 < z < 1.5, and z > 1.5. 
The red solid line represents the local HR04 relation as it would be observed under the biases of our sample. The lines and shaded regions 
in all panels represent the local SMBH mass-bulge mass relation by HR04 and its uncertainty. Our sources are indeed consistent with the 
biased HR04 relation, indicating no significant evolution of the Mbh — M* relation. 



Fig. 9.— Distributions of the specific SMBH and galaxy mass growth rates. The left (right) panel is for the high-(low-)redshift sub¬ 
sample. The specific galaxy growth rate distribution evolves significantly as a function of redshift while the specific SMBH growth rate 
distribution does not change dramatically. Note that 4>* is the probability density function (i.e., f 4>*(cc) dx = 1). 
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o 


Fig. 10.— The “flow patterns” of SMBHs and their host galaxies in the SMBH mass-galaxy total stellar mass plane. Arrows represent 
the direction of the evolution (45 degrees indicates log(sM/sSFR) = 0). Colors indicate the absolute value of the total specific growth rate 

(i.e., yj sM 2 + sSFR 2 ) which spans > 2 orders of magnitude. The dashed line and shaded area represent the local SMBH mass-bulge mass 
relation of HR04 and its uncertainty. The cross indicates the lcr uncertainties of Mbh and M*. T he l ength of the vectors are arbitrary, 
since the “flow patterns” over a long timescale 1 Gyr) depends on AGN duty cycle (see Section |6.3D . 
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Fig. 11.— The same as Figure ITol but for the high-redshift (left panel) and low-redshift (right panel) sub-samples. 
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Fig. 12. — Upper panels: The ratio of the SMBH to the host-galaxy specific mass growth rate vs. the ratio of the SMBH mass to the 
host-galaxy total stellar mass. Left (right) panels are for the high-(low-) redshift sub-samples. The red solid line in each panel represents 
the best fit. The best fit indicates an anti-correlation between and sM/sSFR is significant (null probability p < 10 —6 ). The red 

vertical dashed line in each panel represents the local mass ratio from HR04. The red horizontal dashed line in each panel is an indicator of 
sM/sSFR = 1, i.e., the specific SMBH accretion rate and the specific SFR are equal. Note that there are more sources above the horizontal 
line. This asymmetric distribution indicates that, if the growth rates of our sources persist over some length of time, AGN duty cycle 
is less than unity (i.e., the lifetime of an active SMBH is smaller than the lifetime of star formation), otherwise the SMBH will become 
over-massive as we evolve our sources. Lower panels: The ratio of the SMBH accretion rate to SFR vs. the ratio of the SMBH mass to the 
host-galaxy total stellar mass. Again, left (right) panels are for the high-(low-) redshift sub-samples. A linear fit indicates that M /SFR is 
a constant (with scatter) over M-br/M*. 













































































































































Mbh — M* Evolution and AGN Duty Cycle 
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Fig. 13.— The PDF of tsf and AGN duty cycle (which measures the fraction of BLAGNs among star-forming galaxies, see Section [6.3D 
that would result in the largest fraction of our sources lying on the local SMBH mass-bulge mass relation. Left (right) panels are for 
the high-(low-) redshift sub-samples. The contour in each panel indicates the la (i.e., 68.3%) uncertainty. The color map in each panel 
indicates the probability density. 
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TABLE 1 

Summary of Sample Selection 


Selection Criteria COSMOS CDF-S 

Total z < 2.4 X-ray BLAGNs 221 30 

X-ray BLAGNs with optical/UV and Herschel photometry 85 11 

Final sample: X-ray BLAGNs with optical/UV photometry, robust IR luminosity, good spectra 62 7 



Mbh — M* Evolution and AGN Duty Cycle 

TABLE 2 

Summary of Properties of AGNs and Host Galaxies 


ID 

RA 

deg 

DEC 

deg 

z 

SF a 

Field 

log M B h 
Mq 

log M * 
Mq 

log L b ol b 
erg s —1 

log Li R d 
erg s —1 

A (log Lm) 

erg s —1 

1 

149.83875 

2.67511 

0.26 

S 

COSMOS 

7.07 

10.88 

44.51 

44.93 

0.06 

2 

150.15837 

2.13961 

1.83 

Z 

COSMOS 

8.99 

11.75 

46.08 

46.74 

0.07 

3 

150.34592 

2.14758 

1.26 

z 

COSMOS 

8.28 

11.23 

45.38° 

46.17 

0.08 

4 

150.0425 

2.62917 

1.57 

z 

COSMOS 

8.68 

11.90 

46.16 

46.44 

0.05 

5 

149.69879 

2.44122 

1.52 

I 

COSMOS 

8.74 

10.52 

45.77 

46.17 

0.07 

6 

149.881 

2.45083 

1.32 

z 

COSMOS 

8.01 

11.05 

45.28 

46.11 

0.02 

7 

149.99158 

1.72428 

1.62 

s 

COSMOS 

8.19 

11.55 

46.02 

46.23 

0.04 

8 

150.30954 

2.39914 

1.80 

s 

COSMOS 

8.81 

11.49 

46.46 

46.35 

0.04 

9 

150.62525 

1.80289 

0.63 

I 

COSMOS 

7.43 

10.79 

45.90 

45.67 

0.09 

10 

150.05871 

2.47739 

1.26 

z 

COSMOS 

8.58 

11.52 

45.66 

46.03 

0.03 

11 

149.70587 

2.41975 

1.12 

z 

COSMOS 

9.36 

11.37 

45.77 

45.92 

0.03 

12 

149.94171 

2.79544 

1.07 

s 

COSMOS 

8.67 

11.38 

46.19 

45.90 

0.07 

13 

149.58283 

2.48433 

0.34 

s 

COSMOS 

8.15 

11.25 

45.01 

44.76 

0.01 

14 

150.12508 

2.86175 

1.58 

I 

COSMOS 

8.30 

11.80 

46.08° 

45.90 

0.06 

15 

150.32737 

2.46094 

1.05 

I 

COSMOS 

8.70 

11.29 

45.73 

45.75 

0.10 

16 

149.59246 

1.75675 

1.96 

s 

COSMOS 

8.90 

11.23 

47.15 

46.66 

0.26 

17 

150.14708 

2.71747 

1.18 

z 

COSMOS 

8.13 

10.79 

45.05 

45.63 

0.12 

18 

150.19558 

2.00442 

1.92 

z 

COSMOS 

9.02 

11.14 

46.20 

46.22 

0.11 

19 

150.63387 

2.59369 

0.66 

s 

COSMOS 

8.01 

10.12 

45.59 

45.42 

0.06 

20 

150.49567 

2.41256 

1.37 

z 

COSMOS 

8.43 

11.19 

45.37 

45.64 

0.16 

21 

150.57633 

2.18142 

0.55 

I 

COSMOS 

8.29 

11.06 

44.48 

44.93 

0.00 

22 

150.64521 

2.71481 

0.20 

s 

COSMOS 

6.55 

10.63 

43.81 

44.07 

0.04 

23 

150.19504 

1.79383 

1.87 

I 

COSMOS 

8.74 

11.57 

45.69 

46.09 

0.10 

24 

149.41992 

2.03553 

1.48 

I 

COSMOS 

9.12 

11.49 

45.74 

45.84 

0.10 

25 

150.58187 

2.28769 

1.34 

z 

COSMOS 

7.99 

10.83 

45.72 

45.90 

0.09 

26 

149.47792 

2.64247 

1.6 

s 

COSMOS 

8.64 

11.77 

46.26 

45.94 

0.22 

27 

149.86558 

2.00306 

1.25 

z 

COSMOS 

8.11 

10.41 

45.19 

45.47 

0.10 

28 

149.64183 

2.74089 

1.89 

s 

COSMOS 

8.72 

10.18 

46.45 

46.03 

0.24 

29 

150.658 

2.7835 

0.21 

I 

COSMOS 

6.95 

10.66 

44.85 

44.29 

0.04 

30 

150.55487 

2.641 

1.14 

z 

COSMOS 

8.51 

10.76 

45.33 

45.14 

0.35 

31 

149.58521 

2.05111 

1.35 

z 

COSMOS 

8.80 

11.28 

45.93 

45.77 

0.29 

32 

150.19975 

2.19089 

1.51 

z 

COSMOS 

8.89 

11.18 

45.56 

45.57 

0.16 

33 

150.11929 

2.85353 

0.77 

I 

COSMOS 

9.09 

11.12 

45.60 

45.12 

0.07 

34 

150.31367 

2.80361 

1.46 

I 

COSMOS 

9.06 

10.80 

45.69 

45.58 

0.25 

35 

149.91692 

2.38522 

1.13 

I 

COSMOS 

8.59 

11.42 

45.68 

45.41 

0.22 

36 

149.78971 

2.32125 

0.38 

I 

COSMOS 

7.39 

10.55 

44.63 

44.14 

0.01 

37 

150.6355 

1.66917 

1.79 

I 

COSMOS 

8.78 

10.33 

45.74 

45.31 

0.62 

38 

150.30004 

2.70936 

0.73 

I 

COSMOS 

7.61 

9.33 

44.50° 

44.77 

0.09 

39 

150.13954 

1.6365 

0.52 

z 

COSMOS 

7.49 

10.30 

44.53 

43.93 

1.41 

40 

150.23083 

2.57817 

1.40 

z 

COSMOS 

8.87 

11.35 

45.87 

46.02 

0.20 

41 

150.80187 

2.00061 

1.78 

I 

COSMOS 

9.21 

10.50 

45.47° 

46.06 

0.07 

42 

150.62771 

2.741 

0.82 

I 

COSMOS 

7.92 

10.50 

45.43 

45.12 

0.10 

43 

150.384 

1.57247 

1.36 

s 

COSMOS 

8.17 

11.09 

45.92 

45.85 

0.11 

44 

150.25267 

2.48642 

2.04 

I 

COSMOS 

9.06 

11.39 

45.89 

46.16 

0.13 

45 

150.24517 

1.90008 

1.56 

s 

COSMOS 

8.66 

11.18 

45.97 

45.97 

0.04 

46 

150.19467 

2.06792 

0.55 

I 

COSMOS 

7.29 

11.17 

44.93 

44.77 

0.08 

47 

150.13908 

1.877 

0.83 

I 

COSMOS 

8.12 

11.07 

45.18 

45.09 

0.01 

48 

150.006 

2.81242 

0.77 

s 

COSMOS 

8.38 

10.77 

45.59 

45.10 

0.10 

49 

149.86796 

2.35192 

0.35 

I 

COSMOS 

6.84 

10.63 

44.09 

44.51 

0.02 

50 

149.62429 

2.18067 

1.19 

I 

COSMOS 

8.43 

10.99 

45.55 

45.57 

0.13 

51 

149.4795 

2.80183 

1.11 

s 

COSMOS 

8.68 

11.29 

45.97 

45.37 

0.19 

52 

149.868 

2.33075 

1.49 

I 

COSMOS 

8.67 

11.18 

45.67 

46.10 

0.07 

53 

150.55046 

1.709 

0.37 

I 

COSMOS 

7.35 

10.44 

44.49 

44.34 

0.08 

54 

149.6692 

2.07406 

0.34 

z 

COSMOS 

8.31 

10.67 

44.42 

44.17 

0.10 

55 

150.10521 

1.98117 

0.37 

s 

COSMOS 

8.12 

8.75 

45.23 

44.18 

0.13 

56 

150.05379 

2.58967 

0.7 

s 

COSMOS 

8.33 

11.21 

45.78 

44.88 

0.11 

57 

150.1237 

2.35825 

0.73 

z 

COSMOS 

7.21 

10.62 

44.83 

44.91 

0.11 

58 

150.68317 

2.57461 

0.38 

z 

COSMOS 

7.75 

10.31 

44.78 

44.07 

0.12 

59 

150.2082 

1.87536 

1.16 

I 

COSMOS 

8.83 

10.74 

45.36 

45.31 

0.14 

60 

149.89484 

2.17444 

1.32 

z 

COSMOS 

8.57 

10.90 

45.19 

45.75 

0.09 

61 

149.76346 

2.33411 

1.13 

z 

COSMOS 

8.08 

10.94 

45.28 

45.37 

0.13 

62 

149.66362 

2.08519 

1.22 

z 

COSMOS 

8.19 

10.80 

45.54 

45.48 

0.13 

63 

53.22033 

-27.85556 

1.23 

S04 

CDF-S 

8.25 

10.74 

45.06 

44.96 

0.07 

64 

53.11037 

-27.67658 

1.03 

S04 

CDF-S 

7.80 

10.64 

45.77 

45.55 

0.02 

65 

53.10483 

-27.70525 

1.62 

S04 

CDF-S 

8.74 

11.00 

45.64 

46.15 

0.04 

66 

53.07146 

-27.71761 

0.57 

S04 

CDF-S 

7.65 

10.35 

44.72 

45.33 

0.14 

67 

53.0455 

-27.73756 

1.615 

S04 

CDF-S 

8.84 

11.06 

45.58 

45.58 

0.01 

68 

53.25642 

-27.76183 

0.62 

S04 

CDF-S 

7.96 

10.82 

44.72 

44.85 

0.03 

69 

53.12525 

-27.75656 

0.96 

S04 

CDF-S 

7.43 

10.67 

44.75 

45.59 

0.02 
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TABLE 2 — Continued 


ID 

RA 

DEC 

z 

SF a 

Field 

log M B h 

log M» 

log L Bo l b 

log ii R d 

A(logL IR ) e 


deg 

deg 




Mq 

Mq 

erg s —1 

erg s —1 

erg s 1 


Note. — 

a Spectrum source flag: “S”, “I” and “Z” indicate the spectrum and redshift are from the SDSS archive, the COSMOS Magellan/IMACS campaign 
dTrump et al.ll2009al) and the zCOSMOS VLT/VIMOS campaign (|Lillv et al.ll2007l) . respectively; “S04” means the spectrum and redshift are from 
iSzokolv et al.l (j2004T) . 

^ M can be calculated from LboI using Equation [3] 

C Objects that are detected in the soft X-ray band but not in the hard X-ray band. 

^ SFR can be calculated from Lir using Equation [6] 

e The uncertainty of the IR luminosity reported here is only relevant to the uncertainty of the FIR SED calibration. To fully account for the 
uncertainty of the IR luminosity, one should also add the intrinsic scatter of the SED template in quadrature (see Section l3.5l) . 












